UNCLASSIFIED 


402307 


DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  In  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 

Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  Implication  or  other¬ 
wise  as  In  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


402307 


( 


o 

oo 

<tt 

o 


CURRENT  TRANSITION  GAS-METAL 


Ilf.  STABILITY  OF  GLOW  AND  ARC  DISCHARGES 


G.  Ecker,  W.  Kroll  and  O.Zfiller 


/  ftp  •  ' 

^Technical  Report  FTR  III  ) 


( 


lostitut  fur  Theoretisdie  Physik  der  Universitdt  Bonn 
'Bonn/Germany 


4 


4 


&> 


CURRENT  TRANSITION  GAS-  METAL. 


III.  STABILITY  OF  GLOW  -  AND  ARC-DISCHARGES, 


* 


Ecker,  W. 


KrOll  and  O.  Zdller. 


Technical  Report  FTR  III 


(P  "Cf/>  •  • 


The  reaearch  reported  in  this  document  haa  been  made 
poaaible  through  the  aupport  and  aponaorahip  of  the  US 
Department  of  Army,  through  ite  European  Reaearch 
Office. 


DA  -  91  -  591  -  EUC  -  1764 


) 


3 


ft 


< 


V 


CONTENTS. 


Page 


LIST  OF  SYMBOLS  5 

I  INTRODUCTION  7 

II  THE  CONCEPT  OF  STABILITY  9 

III  MODE  ANALYSIS  OF  THE  STABILITY  PROBLEM  10 

IV  VARIATIONAL  PRINCIPLE  ANALYSIS  OF  THE 

STABILITY  PROBLEM  18 

V  STABILITY  CONSIDERATIONS  IN  PLASMA  PHYSICS  20 

VI  BASIC  EQUATIONS,  BOUNDARY  CONDITIONS 

AND  APPROXIMATIONS  22 

A  THE  GLOW  APPROXIMATION  24 

B  THE  ARC  APPROXIMATION  25 

C  LIMITATION  OF  GLOW-  AND  ARC 

APPROXIMATION  26 

D  INEFFECTIVE  BOUNDARY  CONDITIONS  27 

VII  STATIONARY  MODES  OF  THE  GLOW-AND 

ARC  DISCHARGE  28 

A  NORMAL  GLOW  COLUMN  (SCHOTTKY  COLUMN)  28 
B  THERMAL  GLOW  COLUMN  31 

C  GLOW  COLUMN  WITH  NEGATIVE  IONS  39 

D  WALL  DETACHED  GLOW  COLUMN  46 

E  WALL  STABILIZED  ARC  COLUMN  47 

F  WALL  DETACHED  ARC  COLUMN  54 

G  CATHODIC  PART  OF  THE  ARC  DISCHARGE  61 

H  ANODIC  PART  OF  THE  ARC  DISCHARGE  68 


I  CATHODIC  PART  OF  THE  GLOW  DISCHARGE  69 

J  ANODIC  PART  OF  THE  GLOW  DISCHARGE  77 

THERMAL  INSTABILITY  OF  THE  GLOW  COLUMN  80 

INSTABILITY  OF  THE  CATHODIC  PART  OF  THE 
GLOW  DISCHARGE  91 

STABILITY  PRINCIPLES  IN  GASEOUS  ELECTRONICS  103 

A  THE  PRINCIPLE  OF  MINIMUM  ENTROPY 

PRODUCTION  103 

B  STEENBECK'S  MINIMUM  PRINCIPLE  106 

MONOGRAPHS  109 


LITERATURE 


109 


5 


LIST  OF  SYMBOLS. 


This  list  contains  only  those  symbols  which  are  used  consistently 
throughout  the  whole  report.  Symbols  which  have  only  special 
meanings  for  certain  sections  will  be  appropriately  defined  in 
the  text. 
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magnetic  field 

diffusion  coefficient 

electric  field 

unit  charge 

distribution  function 

Planck's  constant 

saturation  current 

identity  tensor 

Bessel  functions 

Neumann  functions 

electric  current  density 

Boltzmann's  constant 

mass  of  the  charge  carriers 

particle  density 

pressure  tensor 

scalar  pressure 

radius  of  the  discharge  vessel 

space  vector 

radial  distance 

temperature 

time 

ionization  potential 
electric  potential 
drift  velocity 

ionization  coefficient  per  particle  and  per  sec. 
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dt  ionization  coefficient  per  particle  and  per  cm 
Jf  second  Townsend  coefficient 
36  specific  heat  capacity 
X.  mean  free  path  of  the  electrons 
/JL  mobility 
j’gmass  density 
(S'  recombination  coefficient 
eigenvalue 

f  particle  current  dens-ty 
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I.  INTRODUCTION. 


The  object  of  gaseoue  electronic*  ia  a  aystem  which  conaiata  of  varioua 
particle  componenta.  Some  or  all  of  theae  component*  may  be  aingly  or 
multiply  charged  iona.  It  ia  aaaumed  that  there  ia  no  coordination  of  the 
aame  neighbouring  particle*  for  timea  essentially  longer  than  the  life¬ 
time  of  the  aystem  in  a  microstate  of  the  f  -space. 

Obviously  there  is  an  enormous  variety  of  systems  which  fulfil  these 
conditions.  The  quantities  of  these  systems  depend  strongly  on  the  ex¬ 
ternal  parameters  and  the  boundary  conditions  which  are  imposed  by 
the  experimental  setup.  Nevertheless  the  principal  properties  >~f  the 
problem  of  stability  can  be  discussed  in  general  form. 

Let  us  introduce  a  set  of  parameters  Xp  which  we  will  call  the 
EXTERNAL  PARAMETERS.  In  this  group  we  will  take  all  the  quan¬ 
tities  which  define  the  experimental  setup  and  which  can  be  disposed 
of  by  external  manipulation.  Such  parameters  are  for  example  the  geo¬ 
metry  of  the  device,  the  material  and  temperature  of  the  walls,  the  po¬ 
tential  applied,  the  external  magnetic  field  etc.  Let  us  further  intro¬ 
duce  a  set  of  variables  Yy  which  we  call  the  INTERNAL  VARIABLES. 
These  internal  variables  cannot  be  disposed  of  from  outside  but  for  a  gi¬ 
ven  set  of  external  parameters  X9  they  are  determined  by  the  interior 
mechanism  of  the  system.  Of  course  the  number  of  these  variables  intro¬ 
duced  must  be  sufficient  to  allow  a  complete  description  of  the  mechanism 
within  the  range  of  accuracy  required.  Examples  for  such  internal  vari¬ 
ables  are  e.  g.  the  particle  number  densities  n  ,  the  particle  velocities 
v  ,  the  temperatures  of  the  particle  components  T  ,  the  local  elec¬ 
tric  field  E  ,  etc. 

The  time  and  local  dependences  of  the  internal  variables  are  governed 
by  a  system  of  differential  equations  prescribed  by  the  relevant  physi- 
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cal  lawa.  The  external  parameters  enter  these  equations  and  determine 
the  boundary  conditions. 

For  the  STATIONARY  STATES  this  general  system  of  equations  sim. 
plifies  due  to  the  relation  B/M-o.  In  the  following  we  call  such  sta¬ 
tionary  states  "S -MODES". 

In  general  the  possibility  cannot  be  excluded  that  we  have  more  than 
one  s-mode  for  one  given  set  of  external  parameters.  There  are  two 
reasons.  First,  our  system  with  the  boundary  conditions  can  have  more 
than  one  solution.  Secondly,  one  of  the  boundary  conditions  may  be  ren¬ 
dered  ineffective  due  to  the  statistical  fluctuations  always  present  in  a 
real  physical  system. 

As  an  example  we  may  refer  to  the  S  chottky  solutions  for  different 
currents  belonging  to  the  external  parameters  of  the  longitudinal  field, 
tube  radius,  gas  kind  and  pressure.  Interesting  examples  are  also  pro¬ 
vided  by  Kenty's  diffuse  and  filamentary  discharge  (Fig.  15  )  and  by 
the  various  arc  spot  modes  (Fig.  16). 

Let  us  form  a  multi-dimensional  space  of  the  parameters  Xy  .  Then 
let  us  introduce  a  function  U  (Xy  )  which  can  take  only  integer  va¬ 
lues  (U=0,  1,2....).  The  value  of  this  function  determines  the  number 
of  possible  s -modes  for  the  set  of  external  parameters. 

Independent  of  the  value  of  the  U  function  it  is  quite  possible  that 
the  experiment  will  not  produce  any  s-mode.  The  reason  would  be  that 
all  modes  are  unstable.  That  means  they  are  unable  to  sustain  them¬ 
selves  in  the  stationary  state  against  small  perturbations,  which  always 
occur  in  any  kind  of  a  physical  system.  Here  is  the  point  where  the  pro¬ 
blem  of  stability  comes  into  consideration  and  we  see  immediately  that 
it  is  of  excellent  importance. 
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II.  THE  CONCEPT  OF  STABILITY. 

We  call  a  s -mode  stable  if  small  deviations  from  the  Yy  state  are 
damped.  We  call  it  unstable  if  the  small  deviations  from  the  state  grow 
in  amplitude  in  such  a  way  that  the  system  progressively  departs  from 
the  initial  state  and  never  returns  to  it. 

This  general  definition  is  not  sufficient  for  detailed  discussions  be¬ 
cause  the  following  questions  are  not  taken  into  account. 

1.  The  problem  of  stability  can  depend  on  the  amplitude  of  the  distur¬ 
bance.  A  system  may  be  stable  against  small  deviations  but  un¬ 
stable  for  slightly  larger  deviations. 

2.  The  problem  of  stability  can  depend  on  the  kind  of  disturbance  con¬ 
sidered.  A  certain  s-mode  may  be  stable  with  respect  to  disturban¬ 
ces  of  a  certain  kind  but  unstable  with  respect  to  different  distur¬ 
bances. 

3.  If  a  particular  s-mode  is  stable  in  one  part  of  the  Xy  -space  and 
unstable  in  the  other  part  then  there  must  be  a  region  of  the  para¬ 
meters  Xy  in  which  it  is  neither  stable  nor  unstable  but  neutral. 
The  situation  in  this  "neutral  state"  shows  two  different  forms  de¬ 
pending  on  whether  or  not  the  imaginary  part  of  the  coefficient  de¬ 
termining  the  time  dependence  is  zero.  If  the  imaginary  part  is  zero 
then  in  the  neutral  state  the  eystem  passes  from  an  aperiodically 
growing  to  an  aperiodically  damped  state.  If  the  imaginary  part  is' 
finite  then  the  system  passes  from  an  oscillatory  growing  to  an  os¬ 
cillatory  decaying  state. 

We  base  our  discussions  on  the  following  definitions. 

ABSOLUTE  STABILITY:  A  given  system  in  a  certain  stationary  state 
shows  absolute  stability  if  it  returns  to  the  stationary  state  after  any 
kind  and  any  magnitude  of  parameter  disturbance. 


RELATIVE  STABILITY:  A  given  ayetem  in  a  certain  stationary  slate 
shows  relative  stability  if  it  returns  to  this  stationary  state  after  any 
kind  of  infinitesimal  parameter  disturbance. 

TYPE-A-STABILITY:  A  given  system  in  a  certain  stationary  state 
shows  relative  or  absolute  stability  of  type  A  if  it  returns  to  this  sta¬ 
tionary  state  after  an  infinitesimal  respectively  arbitrary  disturbance 
of  parameter  A. 

MARGINAL  INSTABILITY:  If  all  initial  states  are  classified  as  stable 
or  unstable  according  to  the  criteria  stated  above  then  in  the  P  -space 
the  locus  which  separates  the  two  classes  of  states  defines  the  state  of 
marginal  or  neutral  stability. 

SECONDARY  FLOW  INSTABILITY:  If  in  the  state  of  marginal  instabi¬ 
lity  the  imaginary  part  of  the  coefficient  determining  the  time  dependence 
is  zero  then  we  say  that  the  principle  of  the  exchange  of  stabilities  is 
valid  or  that  we  have  an  instability  of  secondary  flow.  In  this  case  the 
restoring  forces  are  to  small  to  sustain  the  s-mode. 

OVERSTABILITY :  If  at  the  onset  of  the  instability  the  imaginary  part 
of  the  coefficient  determining  the  time  dependence  is  finite  then  the  os¬ 
cillatory  motion  prevails  and  we  have  a  state  of  over  stability.  In  this 
case  the  restoring  forces  are  so  large  that  they  overshoot  the  stationary 
state  of  the  s-mode  and  cause  the  increasing  oscillations. 


III.  MODE  ANALYSIS  OF  THE  STABILITY  PROBLEM. 


In  this  investigation  we  will  not  treat  the  problem  of  "absolute  stabi¬ 
lity". 

The  theory  of  relative  stability  has  the  decisive  advantage  to  allow 
linearization  of  the  equations. 

Our  preceding  definition  of  "relative  stability"  suggests  the  following 
way  to  investigate  the  stability  question: 
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We  consider  small  perturbations  of  the  internal  variables  and  study 
<  the  time  dependence  of  these  deviations.  Here  two  difficulties  occur. 

The  one  is  that  our  system  is  stable  only  if  it  is  stable  with  respect 
to  all  possible  initial  increments.  The  second  is  that  the  time  depen¬ 
dence  not  only  changes  the  amplitude  of  a  given  initial  disturbance  but 
also  its  form.  It  is  possible  that  with  the  change  in  the  form  an  ini¬ 
tially  growing  amplitude  might  gradually  die  out  again.  Such  a  state 
we  would  not  call  unstable. 

To  overcome  the  difficulties  described  above  we  consider  the  eigen- 
solutions  of  our  perturbation  (index  (1)  )  which  are  characterized  by 

the  time  dependence 

pl1\r,t)  •  Fl1>(r)eu*t  (»> 

Introducing  equation  (1)  into  our  system  we  find  a  set  of  equations  which 
defines  together  with  the  boundary  conditions  an  infinite  number  of  ei- 
gensolutions  and  eigenvalues  at  .  These  eigenfunctions  provide  a 
complete  orthogonal  system  of  functions.  Any  initial  disturbance  may 
be  developed  after  these  functions.  Consequently  the  question  of  stabi¬ 
lity  is  reduced  to  the  question  whether  anyone  of  the  eigenvalues 
has  a  positive  real  part. 

Therewith  our  task  is  to  determine  the  eigenvalues. 

Only  in  special  cases  the  eigenfunctions  can  be  determined  in  closed 
form  and  we  have  directly  a  dispersion  relation 

WM  •  0  (Z) 


defining  the  eigenvalues  ou 


In  general  normal  mode  analysis  is  applied.  In  this  case  the  eigen- 
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function  is  developed  into  a  series  of  appropriate  orthonormal  functions, 
(normal  modes). 

(3)  FM(r)  •  Zakfk(r) 

The  terms  of  the  differential  equations  can  still  depend  on  the  coordi¬ 
nates  via  the  imperturbed  stationary  solution.  Therefore  we  encounter 
products  of  the  form 

(4)  g<0>(r)ZakfJF ) 

Now  we  use  for  this  whole  term  a  mode  representation  in  the  form 

< 5)  g  <•>(?)  ZaJM(r)-l  b,  f,  Ir) 

k  j  J  J 

where  the  coefficients  bj  are  related  to  the  coefficients  a^  by 

(6)  bj  -  Z  [ JgM(r)  fk  (r~)  fj  (?)  w  (r)  dz ]ak 

We  introduce  this  into  our  equations  and  collect  all  terms  belonging 
to  the  same  normal  mode.  Since  our  relations  must  be  fulfilled  for  any 
value  of  the  coordinates,  the  sum  over  the  coefficients  of  each  mode 
must  be  identical  zero  for  all  equations.  This  provides  an  infinite  set 
of  homogeneous  equations.  We  have  as  many  equations  as  coefficients 
a^  and  therefore  a  solution  exists  only  if  the  secular  determinant  of 
the  coefficients  disappears.  This  condition  leads  to  the  characteristic 
dispersion  equation  for  the  eigenvalues  (jj 

Introducing  these  eigenvalues  into  the  equations  for  a^  allows  to 
calculate  these  coefficients  and  with  that  the  eigenfunctions. 

If  we  have  already  the  eigenvalues  u>  from  equation  (2)  or  from  the 
secular  determinent  then  we  can  decide  the  question  of  stability  by 
studying  the  sign  of  the  real  parts  of  u> 
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However  if  we  are  merely  interested  in  stability  we  do  not  really  have 
to  solve  the  dispersion  relations.  The  question  whether  all  the  roots 
lie  in  the  negative  Gaussian  plane  can  be  answered  by  the  following 
criteria. 

Following  Nyqui  st  let  us  assume  that  our  general  dispersion  rela¬ 
tion  (2)  describes  a  conformal  mapping  of  the  Ui  -plane  in  the 
W-plane.  Then  from  Cauchy"  s  theorem  we  have 

W(Ui)  dta  •  27  i  N  (7) 

X 

Where  N  is  the  number  of  zeros  of  W  within  the  contour  of  the 
integral.  Using  the  substitution  W  =  R.  exp  (i$  )  it  follows  then 

Af  -  27  iN  (8) 

Consequently  the  question 
of  stability  may  be  decided 
by  plotting  in  the  W-plane 
the  curve  corresponding  to 
the  integral  contour  in  the 
U)  -plane  (See  Fig.  1). 

From  the  number  of  en- 
circlings  of  this  curve  we 
may  conclude  to  the  num¬ 
ber  of  roots  lying  within 
the  contour  of  our  integral 
in  the  (0  -plane.  This 
provides  us  with  the  cri¬ 
terion  for  stability. 

In  the  case  of  the  normal  mode  analysis  the  dispersion  relation  is  in 
general  represented  by  an  infinite  power  series. 


n$.  it  aiMMUMi  mui  s  uw  mnH' 


H 


It  can  happen  that  this  series  reduces  to  a  polynomial  or  can  be  repre¬ 
sented  as  a  product  of  polynomials.  Then  we  have  simpler  criteria  to 
decide  the  question  of  stability. 

A  necessary  condition  for  all  the  roots  having  negative  real  parts  -that 
means  a  necessary  condition  for  stability  -  is  that  all  the  coefficients 
of  the  polynomial  have  the  same  sign. 

If  the  polynomial  considered  is  of  the  quadratic  type  then  the  criterion 
of  equal  sign  of  all  coefficients  is  also  sufficient.  For  higher  order  po¬ 
lynomials  the  sufficient  criteria  are  more  complicated. 

The  best  known  necessary  and  sufficient  criteria  to  decide  whether  all 
roots  have  negative  real  parts  are  those  ofRouth  and  Hurwitz. 

According  to  the  criterion  of  Hurwitz  we  have  stability  provided  that 
all  elements  c  and  all  "Hauptabschnittdeterminanten"  D  of  the 
following  matrix 

c,  Co  o  o  •  •  •  o 

c §  c*  cf  cQ  •  o 

CS  C4  Cjf-  •  o  > 

•  • 

•  • 

OOOO  cn 

d 

have  a  positive  sign.  The  elements  of  the  matrix  correspond  to  the  co¬ 
efficients  of  the  secular  equation.  The  "Hauptabschnittsdeterminenten" 
are 

»•’  O’-*. 


(9)  j 


•  •  • 


m  Cn  Dfi-t 
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The  Routh  criterion  can  be  understood  from  the  following  scheme. 


Cn.rc„.s  Cn-3  •  •  '(-c„)/cn.,  ^ 

an-4  •  ’  • 

bn-Mbn-3  *n-f  ‘  ’  ’  Y* 

•  •  »  * 

•  •  •  • 

•  •  •  .  .  •  • 

In  this  scheme  the  third  equation  is  produced  from  the  first  two  ones 
by  multiplying  the  second  with  the  factor  indicated  on  the  right  hand 


side  and  adding  it  to  the 
first  equation.  The  result 
is  given  in  the  third  equa¬ 
tion  where  we  have  shifted 
every  coefficient  by  one 
place  to  the  left  hand  side. 
The  fourth  equation  is 
produced  in  a  similar  way 
by  the  multiplication  of  the 
third  equation  and  addition 
to  the  second  one  etc.  If 


all  the  roots  are  supposed  .i..«k t.r •». »or,.»Ur. 

vsnkrllerium", 

to  have  negative  real  parts, 

then  all  the  numbers  of  the  first  column  must  be  positive.  That  means 
proceeding  from  one  equation  to  the  next  we  may  stop  our  calculation 
at  the  very  moment  if  zero  or  a  negative  quantity  appears,  because  then 
we  know  that  the  problem  under  consideration  exhibits  instability. 


The  "ORTSKURVENKRITERIUM"  considers  the  mapping  of  the  com¬ 
plex  U)  -plane  into  the  W-plane  as  prescribed  by  the  dispersion  re¬ 
lation  W(  U)  ).  The  imaginary  axis  of  the  U)  -plane  (  u>  =  iy)  pro¬ 
duces  a  certain  curve  in  the  W-plane  (See  Fig.  2). 
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IW(uf)  •  C0  *  CrC0  *  CfUf* *  •  •  •  t  Cff  Oin 

-  Ca  ’CMy  +  c4y*-  +  --  +  i(c,y-c4y*+csy*-  + 
•  X(y)  *  i  /(/) 

-Re'* 


If  we  follow  the  imaginary  axis  of  the  u>  -plane  from  -s*to  +  oo 
then  again  a  certain  path  in  the  W-plane  the  socalled  "Ortskurve"  is 
travelled  trough.  It  is  sufficient  to  consider  only  the  part  /<+  00 
since  we  have  the  relations  X(-  y)  =  X  (y)  and  Y  (-y)  =  -  Y(y). 

The  criterion  reads  now:  Equation  (12)  has  roots  with  negative  real 
parts  only  if  the  "Ortskurve"  in  the  W-plane  surrounds  the  origin 
in  such  away  that  the  angle  $  of  the  pointer  in  the  W-plane  covers 
the  range  4-0 . nT/i  if  y  goes  from  zero  to  infinity. 

If  there  are  roots  with  positive  real  parts  and  the  angle  covered  may 
be  written  in  the  form 

(13)  A  $  -  (q.p)I.  m  tJT 

then  the  number  of  roots  with  positive  real  part  is  given  by 

O4)  p  •  Zj* 

and  number  of  roots  with  negative  real  part  is  given  by 

(15)  q  -  SZl 

To  evaluate  the  "Ortskurvenkriterium"  it  is  not  really  necessary  to 
plot  the  path  in  the  W-plane.  It  is  already  sufficient  to  show  that  the 
following  is  correct: 


■  £  •->. 
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"LUCKENKRITERIUM":  The  equation  W  (  U)  )  =  0  has  roots  with 
negative  real  parts  only  if  the  equations 

X (s)  -o.  ±  Ym  Y( 4)  •  o  (16) 

2 

with  s  =  y  have  positive  roots  s^  respectively  s^  only, if  the  zeros 
follow  the  sequence 

^  ^  ^  Sjf  K  S j  <  •  (17) 

It  is  not  even  necessary  to  solve  for  these  roots  since  it  is  sufficient 
to  show  that  for  c^  >  0  the  equation 


Yfs)  •  o 


(18) 


has  positive  real  parts  only 

O  <sr  <SM  <Sj  "  •  •  (19) 

and  that  the  values  of  the  function  Xj  =  X(sj  )  changes  its  sign  from 
one  to  the  next  root. 


X  f  ^  o  j  Xj(  >  o/  X  j  ^ 


(2o) 


The  normal  mode  analysis  as  described  above  has  been  frequently 
used.  For  a  specific  problem  it  is  a  powerful  instrument.  It  allows  to 
decide  not  only  the  question  of  stability  but  also  implies  the  determi- 

I  ♦ 

nation  of  the  eigenfunctions  and  the  growth  of  the  instability  respective- 
,  ly  the  decay  time. 

’  The  second  procedure  to  tackle  the  problem  of  stability  restricts  it¬ 

self  from  the  outset  to  the  calculation  of  the  stability  criterion.  There- 
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fore  the  mathematical  handling  is  simpler.  The  approach  we  refer 
to  is  the  METHOD  OF  VARIATIONAL  PRINCIPLE. 


IV.  VARIATIONAL  PRINCIPLE  ANALYSIS  OF  THE 
STABILITY  PROBLEM. 

There  are  essentially  two  difficulties  inherent  in  the  normal  mode  ana¬ 
lysis.  After  choosing  an  appropriate  set  of  orthogonal  modes  to  expand 
the  eigenfunctions  the  one  problem  arises  from  the  calculation  of  the 
coefficients  of  the  expansions  (  6  ).  These  coefficients  include  difficult 
integral  expressions  which  cannot  be  evaluated  in  general  form.  The 
other  problem  is  caused  by  the  formidable  difficulties  in  the  evaluation 
of  the  secular  equation. 

A  simplification  can  be  expected  by  renouncing  the  specific  knowledge 
of  the  frequencies  and  growth  rates  and  by  simply  restricting  our 
question  to  the  criterion  for  the  onset  of  instabilities. 

This  is  the  basic  idea  of  the  variational  approach  to  the  stability  problem. 

The  principal  features  of  the  calculations  are  quite  similar  to  those  of 
the  normal  mode  analysis.  We  again  consider  small  deviations  from  an 
s -mode.  This  allows  us  to  linearize  the  relations  and  we  have  a  set  of 
differential  equations  for  the  first  order  perturbation  where  the  coeffi¬ 
cients  may  depend  on  the  qualities  of  the  stationary  state  solution. 

As  in  the  case  of  the  normal  mode  analysis  the  eigenfunctions  of  this 
first  order  problem  are  characterized  by  the  relation  (1). 

We  reduce  our  system  of  differential  equations  to  one  equation 


(21) 


&(yr) m  o 


t 

i 


But 


I 


l 


id 


for  the  dependent  variable  .  We  develop  TjT  after  the  eigen¬ 

functions  ^  in  the  form 

The  applicability  of  the  variational  method  depends  on  the  character 
of  the  operator  X  .  Therefore  specification  of  this  operator  is  required 
for  further  discussion.  As  an  example  we  treat  here  the  case,  where 
equation  (21)  reduces  to 


m  <^k  Yh  (23> 

where  the  operator  does  not  contain  the  time  anymore.  Introducing 

equation  (22)  into  equation  (21)  and  using  equation  (23)  we  find 


2! aktCr  ( y/  )  -  21  o  at  w 

k  "  k  *  H  '  k 


(24) 


We  now  develop  the  adjoint  function  ^  after  the  adjoint  eigenfunctions 

%■ 


(25) 


Multiplication  of  equ.  (24)  with  Tjf  and  the  appropriate  weight  function 
w  and  integration  produces 


Jfk  °j  %  a„*r  dr  ■  fZSj oka>k  fj  yM  w(?)drw 

If  the  operator  X>^  is  self-adjoint  it  follows 


ofcJCr  (y)  w(r)dT 


Z  a  a  cl) 

k  k  k  k 


(27) 
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Remembering  that  the  ere  first  order  quantities  we  aim  to  re¬ 

present  the  right  hand  side  of  equ.  (27)  as  a  variation  of  a  physical 
quantity  Q.  If  this  is  possible  we  arrive  at 

(28)  £a  a  w  ■  &Q. 

k  k  k  k 

Since  all  values  of  a^a^  are  positive,  the  left  hand  side  can  be  ne¬ 
gative  only  if  at  least  one  of  the  eigenvalues  has  a  negative  real  part. 

In  this  simple  example  the  imaginary  parts  must  all  be  zero,  since  the 
physical  quantity  Q  is  real. 

Therefore  if  all  the  eigenvalues  are  real  and  positive,  the  right 

hand  side  of  the  equation  (28)  is  a  positive  definite  form  of  the  devia¬ 
tions  from  the  stationary  state.  Since  it  is  sufficient  if  one  value  be¬ 
comes  negative  we  can  conclude  the  following  stability  criterion: 

The  system  under  consideration  is  Btable  then  and  only  then  if  the  phy¬ 
sical  quantity  Q  is  a  positive  definite  form  with  respect  to  all  pos¬ 
sible  deviations  from  the  stationary  state. 


V.  STABILITY  CONSIDERATIONS  IN  PLASMA 
PHYSICS. 


The  methods  prescribed  in  the  preceding  paragraphs  were  originally 
developed  in  connection  with  the  stability  problems  in  hydrodynamics. 

During  the  recent  years  they  have  been  widely  applied  also  to  problems 
of  plasma  physics.  The  applications  may  be  grouped  according  to  the 
underlying  physical  model. 

The  MHD  approximation  is  based  on  the  essential  assumptions  of 
skalar  pressure,  an  adiabatic  invariant 
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and  infinite  conductivity. 

For  this  model  the  gross  instabilities  of  a  cylindrical  plasma  column 
have  been  calculated  .  Also  the  stability  of  a  plasma  in  contact 
with  a  magnetic  wall  has  been  studied.  Here  normal  mode  analysis 
produces  similar  to  the  R  a  y  1  e  i  gh  Taylor  instabilities  of  hydro¬ 
dynamics,  gravitational,  flute  and  ripple  instabilities 

The  variational  analysis  has  been  applied  to  m^gnetohydrodynamics 
by  Lundquist  and  later  by  B  e  r  n  s  t  e  i  n  ,  Frieman,  Krus- 
kal  and  Kuls  rud  ^  .  They  showed  that  for  a  magnetohydrodynamic 
device  the  function  Q  can  be  identified  with  the  potential  energy. 

That  means  a  magnetohydrodynamic  system  is  in  equilibrium  if  the 
variation  of  the  potential  energy  is  a  positive  definite  quadratic  form, 

-  provided  that  certain  conditions  are  satisfied. 

Chew,  Goldberger  and  Low  ^  have  treated  a  different  appro¬ 
ximation  in  which  the  pressure  is  not  scalar  and  there  are  two  adia¬ 
batic  invariants  for  the  pressure  parallel  and  perpendicular  to  the  mag¬ 
netic  field.  In  the  C.  G.  L.  approximation  the  stability  of  a  plasma  may 

a 

be  also  decided  by  an  energy  variation  7 

Generalizations  of  the  energy  principles  of  the  M.H.  and  the  C.  G.  L. 

approximation  have  been  given  by  K  r  u  s  k  a  1  and  Oberman^  and 

72 

Rosenbluth  and  Rostocker  .  The  energy  change  calculated 

by  these  authors  is  bounded  below  by  the  M.H.  energy  change  and  above 
by  the  C.  G.  L.  energy  change. 

In  a  recent  publication  R  o  s  e  nb  luth  and  Fuerth  study  the 
problem  of  micro  instabilities  in  a  system  of  finite  resistivity. 
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For  the  general  phenomenon  of  micro  inetabilitiea  there  ia  no  general 
phyaical  model  but  varioua  apecial  caaea  have  been  treated  71,  73  . 

For  our  inveatigation  of  the  atability  qualitiea  of  the  arc  and  the  glow 
diacharge  the  reaulta  quoted  above  are  not  applicable,  becauae  the 
underlying  conditiona  are  not  met. 

Since  the  equationa  describing  the  system  are  characteristic  and  de¬ 
cisive  for  the  stability  consideration  one  would  like  to  define  our  pro¬ 
blem  by  stating  the  equations  and  the  approximations. 

Unfortunately  in  the  case  of  the  gas  discharges  there  is  no  unique  set 
of  equations  describing  the  whole  arc  and  glow  region.  Rather  we  have 
in  each  of  these  discharge  types  a  number  of  model  zones  governed  by 
different  relations. 

As  a  consequence  we  will  expect  that  the  various  model  regions  of  the 
discharges  can  exhibit  different  stability  qualities.  Our  discussion  must 
be  specified  with  respect  to  each  model  region. 

Since  all  stability  calculations  require  the  knowledge  of  the  s -modes 
we  at  first  discuss  the  possible  s-modes.  This  presents  at  the  same 
time  the  equations  and  approximations  of  our  considerations. 

VI.  BASIC  EQUATIONS.  BOUNDARY  CONDITIONS. 

AND  APPROXIMATIONS. 

Most  of  our  investigations  will  be  based  on  the  following  equations. 

From  the  Boltzmann  transport  equation  we  derive  by  multiplication 

2  , 

with  m,  mv,  mv  /2  and  integration  over  the  momentum  space  the 
following  transport  equations 


mrn&m0* flMMNM iS&f®# 


57  "  f  v  ("*<!)•  *, 


j:r-£E*r’£*i(r-ri)-*rni.v(Z) 


7-  <WJ*  7  (nv*?)-Z2±i.f.  . 


where  the  approximations 


P  -  nlcTI 


/  f/vjf,  (Vf(Y)ciZ{.rt  dvtW-^Zhit-\)Jix 

have  been  used. 

denotes  the  net  particle  production  per  second  and  unit  volume.  The 
collision  frequencies 

r  (34) 

V  - nj $cdQ 

are  related  to  the  mobilities  and  the  electron-ion  interaction 

parameters  according  to 

v  mSLu  .  i>  .  -0  »  (35) 

^  9  /  i  -f  m.  (  /  f*  ( 


i 


24 


A  THE  GLOW  APPROXIMATION. 

Characteristic  for  the  glow  approximation  is  the  following: 

In  the  particle  conservation  law  for  the  charged  particle  components 
the  destruction  is  dominated  by  diffusion  and  mobility  losses  and 
not  by  volume  recombination.  In  other  words:  The  fact,  that  each 
volume  element  is  an  open  system  with  respect  to  the  charged  par¬ 
ticle  balance  is  decisive  and  characteristic  for  the  glow  approxima¬ 
tion. 

Characte  ristic  for  the  glow  model  is  also  the  boundary  condition 
at  the  wall.  We  have  at  the  sheath  edge  for  the  particle  number  den¬ 
sity  the  condition 

(38)  (0f  grad  n)„  -  (»  ^  )„  '  Xn 

where  is  an  uncertainty  factor  of  magnitude  one,  accounting  for 

the  influence  of  the  geometry  and  transition  region  a  ,  For  very 
small  values  of  the  mean  free  path  equation  (38)  simplifies  to 


Apart  from  these  characteristic  features  one  in  general  also  assumes 
weak  ionization 

(36)  nnkTn  -  p  -  const 

quasineutrality 

n_  •  nf  •  n 


(37) 


.#vr  ■ 
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and  the  validity  of  the  energy  balance* 

T,  •  T„  <«> 

and 

71  •  c  £  X.  HD 

with  C  being  a  constant  typical  for  the  gas. 

B  THE  ARC  APPROXIMATION. 

Characteristic  for  the  arc  approximation  is  the  following: 

In  the  particle  conservation  law  for  the  charged  particle  components 
the  destruction  is  dominated  by  volume  recombination  and  not  by  dif¬ 
fusion  and  mobility  losses.  In  other  words:  The  fact,  that  each  volume 
element  is  a  closed  system  with  respect  to  the  charged  particle  balance 
is  decisive  and  characteristic  for  the  arc  approximation. 

Characteristic  for  the  arc  model  is  also  the  boundary  condition 
at  the  wall 


(Q)w  •  0  (42) 

From  these  conditions  the  local  validity  of  the  Saha  equation 


2p*(kT)k  (2  7m.  >* 

""  - T& - 


(43) 


follows. 
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Apart  from  these  characteristic  conditions  one  in  general  also 
assumes 


(44) 


r* 


r_  •  r* 


r 


and  quasineutrality . 


C  LIMITATION  OF  GLOW  AND  ARC  APPROXIMATION. 

The  glow  and  the  arc  approximation  described  in  the  preceding 
chapters  are  complementary.  We  therefore  need  only  one  limitation. 

According  to-  the  definitions  this  limitation  is  given  by 


(45) 


dir  ( n  j  •  < Y  n  n, 


where  (  G*  )  is  the  recombination  coefficient. 
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D  INEFFECTIVE  BOUNDARY  CONDITIONS. 

The  above  definitions  of  the  glow  and  the  arc  model  are  based  on  the 
assumption  that  the  boundary  conditions  (38)  and  (42)  are  the  effec¬ 
tive  ones. 

Physically  these  boundary  conditions  formulate  the  requirement  of 
particle  current  continuity  and  energy  current  continuity  respectively, 
at  the  sheath  edge.  There  may  be  physical  influences  not  included  in 
our  description  above  which  take  an  effect  on  these  continuity  require¬ 
ments.  These  effects  may  either  be  due  to  additional  phenomena  like 
external  fields,  convection  etc.  or  they  may  be  due  to  the  normal 
statistical  fluctuations  present  in  any  kind  of  real  discharge. 

If  these  influences  become  stronger  than  the  terms  included  in  the 
formulation  of  our  boundary  conditions  (38)  and  (42)  then  they  render 
these  conditions  ineffective. 

The  behaviour  of  our  discharge  is  different  depending  on  whether  the 
ineffectiveness  is  caused  by  an  other  external  influence  or  by  statisti¬ 
cal  fluctuations. 

In  the  first  case  the  additional  phenomenon  determines  the  selection 
of  the  observed  mode  and  we  have  a  different  effective  boundary  con¬ 
dition. 

In  the  second  case  the  statistical  fluctuations  which  outrule  the  boun¬ 
dary  conditions  allow  transitions  between  all  those  modes  which  are 
possible  without  the  wall  boundary  condition.  In  this  case  the  observed 
mode  is  determined  by  a  stability  criterion. 
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VII  STATIONARY  MODES  OF  THE  GLOW-  AND 
ARC  DISCHARGE. 


Without  a  claim  to  completeness  we  underlay  our  investigations  the 
following  stationary  modes  of  the  arc-  and  glow  discharge. 


A)  NORMAL  GLOW  COLUMN  (SCHOTTKY -COLUMN 

B)  THERMAL  GLOW  COLUMN 

C)  GLOW  COLUMN  WITH  NEGATIVE  IONS 

D)  WALL  DETACHED  GLOW  COLUMN 

E)  WALL  STABILIZED  ARC  COLUMN 

F)  WALL  DETACHED  ARC  COLUMN 

G)  CATHODIC  PART  OF  THE  ARC  DISCHARGE 

H)  ANODIC  PART  OF  THE  ARC  DISCHARGE 

I)  CATHODIC  PART  OF  THE  GLOW  DISCHARGE 

J)  ANODIC  PART  OF  THE  GLOW  DISCHARGE 


A.  NORMAL  GLOW  COLUMN  (SCHOTTKY  COLUMN). 


The  normal  glow  column  is  governed  by  the  following  mechanism: 

Electrons  and  ions  move  in  the  axial  direction  under  the  influence  of 
a  longitudinal  electric  field.  At  the  same  time  they  move  towards  the 
wall  under  the  influence  of  a  concentration  gradient  and  a  radial  ambi- 
polar  field.  Particle  production  is  due  to  collisions  of  electrons  with 
neutral  particles.  Particle  loss  is  due  to  wall  recombination  only. 
Volume  recombination  is  neglected.  The  temperature  of  the  electrons 
is  defined  by  the  energy  gain  in  the  electric  field  on  the  one  hand  and  the 
energy  loss  by  coll.sions  with  neutral  particles  on  the  other  hand.  The 
temperature  of  the  ions  is  determined  by  the  energy  gain  from  the  elec¬ 
tric  field  and  the  energy  loss  due  to  collisions  with  the  neutral  particles. 
There  is  no  radiative  loss.  Therefore  the  total  energy  gained  from  the 
l  electric  field  is  carried  to  the  walls  by  heat  conduction  of  the  neutral 

I 
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particles.  The  energy  input  is 
so  small  that  the  neutral  par¬ 
ticle  temperature  can  be  con¬ 
sidered  constant  across  the 
discharge.  Then  -  for  weak 
ionization  -  also  the  electron- 
and  ion  temperature  are  con¬ 
stant.  The  longitudinal  elec¬ 
tric  field  adjusts  itself  in 
such  a  way  that  the  number 
of  particles  produced  in  the 
discharge  covers  the  particle 
loss  to  the  walls. 


PIr.  i>  Schemalical  sketch  of  the  density  ami 
temperature  distribution  in  the  normal 
alow  column. 


Under  these  assumptions  the  glow  approximation  (VI  A)  is  applicable 
throughout  the  whole  discharge  volume. 

The  longitudinal  electric  field  is  constant  due  to 

rot  E  -  0  (46) 

In  equations  (31)  we  use  the  conditions  of  stationarity  a/Bt-o 

,  weak 

ionization  P=nnkTn  and  the  relation 


°t  '  At  Pt/»* 

With  the  condition  of  congruence 


(47) 


(48) 


we  eliminate  all  azimuthal  components  and  the  electric  field  and  arrive 
at 
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(49) 


with 


n 


(50) 


_ l - .  .  o  — i - , 

M.  */t*  1  +/tt/A.  8’  am  I*/**/1-  & 


Introducing  this  into  the  particle  conservation  law  we  have  the  equation 


(51) 


+  oCn 


O 


which  together  with  the  boundary  condition  (39)  defines  an  eigenvalue 
problem  for  the  electron  density  and  the  eigenvalue  efc  .  Since  the 
particle  density  n  can  only  assume  finite  positive  values  the  well 
known  solution  of  this  eigenvalue  problem  is 

(52)  n(r)  -  na  J0  ([*/D9ff]*r) 

with 

(53)  CC  «  (i'U/R)*  Oeff 

Since  the  production  coefficient  (C  depends  on  the  electron  tempe¬ 
rature  according  to  the  relation 

(M)  «C  (T.)  ■  A  (T.)  exp  I-  ) 

and  the  electron  temperature  is  determined  by  the  electric  field 

according  to  equation  (41)  the  eigenvalue  (53)  defines  the  eigenvalue 

s  ... 

}  of  the  electric  field.  Obviously  the  gradient  E  is  independent  of 

f1  * 
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the  current  within  the  frame  of  the  Schottky  theory. 


Schottky  's  theory  underlies  the  assumption  that  the  charged  car¬ 
riers  are  produced  by  a  single  ionizing  collision.  In  rare  gases  where 
the  formation  of  metastable  atoms  is  common  cumulative  ionization  may 
occur.  One  expects  these  processes  to  increase  with  current.  The  cu¬ 
mulative  ionization  causes  a  field  decrease  with  increasing  pressure. 

22  81 

Spenke  and  Fabrikant  ’  included  cumulative  ionization  by 
adding  nonlinear  terms  in  the  equation  (51).  In  an  exact  treatment 
the  additional  particle  components,  like  e.  g,  the  metastables,  should 
be  treated  with  a  separate  particle  balance  equation  (see  e.  g.  Chapter 
VII  C). 

Schottky 's  theory  also  does  not  account  for  volume  recombination. 
This  effect  becomes  important  with  increasing  particle  densities.  Par¬ 
ticularly  in  the  presence  of  electronegative  gases  or  molecular  gases 
volume  recombination  may  be  of  remarkable  influence. 


B  THERMAL  GLOW  COLUMN. 

One  of  the  characteristic  approximations  for  the  Schottky  theory 
of  the  normal  glow  column  is  the  assumption  of  constant  electron,  ion 
and  neutral  gas  temperature  across  the  discharge. 

In  this  chapter  we  will  give  an  approximate  calculation  of  the  column 
in  which  the  thermal  effects  cannot  be  considered  as  a  small  pertur¬ 
bation  of  the  normal  glow  theory. 

For  the  thermal  glow  column  the  whole  model  concept  of  the  normal 
glow  column  holds  with  one  exception:  The  energy  input  is  so  large 
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that  a  variation  of  the  neutral 
gae  temperature  is  caused  which 
cannot  be  considered  as  a  small 
perturbation.  As  a  consequence, 
electron  and  ion  temperature 
show  a  radial  variation. 


Therefore  the  following  equa¬ 
tions  apply. 

The  particle  conservation  law 
for  the  electrons 


rig.  4t  Schematics!  sketch  of  the  density  and 
temperature  diet  ri  but  ten  in  the  thermal 
glow  column. 


L  4  (r4s)  +  -£-n-  O 

r  dr  \r  dr  J  f  Orf 


Under  the  assumption  of  weak  ionization  the  particle  balance  for  the 
neutral  component  simplifies  to 


nn  k  Tn  •  p  •  const 


The  energy  balance  of  the  electrons  reads 


r.  •  c  £  A. 


and  the  energy  balance  of  the  neutral  particle  component  is 


-*-• 4:  (r  $)*•{/>.*/*.) £‘n  *  ° 


»»*..«  *■-«■—  ***-•» 
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The  coefficient  of  the  volume  particle  production  per  sec.  will  be  re¬ 
presented  in  the  form 

oC  -  A  exp  )  (59) 

The  experiments  provide  the  information 

y-  -  A  exp  (-^)  (60) 

Here  it  should  be  remembered  that  the  characteristic  parameters  are 
actually  not  06  /p  and  E/p  but  et/n^  anc*  E/nn  w^ere  n  *8 
the  neutral  particle  density.  Since  all  experiments  are  carried  out  with 
constant  neutral  particle  temperature,  p  reflects  the  dependence  of 
n^  .  In  our  calculations  however  the  temperature  of  the  neutral  gas 
can  cause  a  variation  in  the  neutral  gas  density.  Therefore  we  should 
apply  the  relation 


The  experimental  knowledge  (60)  not  only  provides  the  coefficient  A 
but  also  defines  the  coefficient  C  of  equation  (57)  according 

to 

T  .  ~L  h  c  _  }  e 
T-  kBpT^£-eA-£  (“I 

Now  we  introduce  equations  (56)  and  (57)  into  (58)  and  have  the  following 
relation  for  the  electron  temperature 

)*»(/*,•/*■  )£‘n  *  0 


(63) 
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With  the  abbreviations 


Equations  (55)  and  (63)  are  two  simultaneous  differential  equations  for 
the  electron  temperature  and  the  electron  density  in  the  discharge. 


o 

R 


This  is  an  eigenvalue  problem  with  the  eigenvalue  E 


We  base  our  approximation  on  the 
fact  that  the  production  coefficient 
due  to  the  strong  dependence  on 
the  electron  temperature  shows  a 
shape  similar  to  that  of  a  step 
function.  Consequently  we  repre¬ 
sent  A  by  the  step  function 
shown  in  Fig.  5.  The  three  quan¬ 
tities  ,  «tg  ,  £  ,  characte¬ 

rize  the  step  function,  o Ca  is  the 
production  coefficient  correspon¬ 
ding  to  the  electron  temperature 
in  the  axis.  is  the  production 
coefficient  corresponding  to  the 
wall  temperature  of  the  electrons, 
is  derived  from  the  condition  c6. 


o  ft*  — —  I 


Fig.  ti  Step  approximation  of  the  production 
coefficient. 

The  position  of  the  step  (  f  ) 
/  e6  (  J )  ■  6  to  be 


(66) 


Jo  ¥ )  “  1  ~9Ui  e{/t4  +^.)E*n0  ~oiff 


With  this  step  function  for  s6  it  is  not  difficult  to  find  the  density 
and  temperature  distributions  in  the  interior  and  in  the  outer  part 
independently.  The  result  is 


ni  "  no  Jo  l[*o/09ff  ]  *r) 


(67) 


"  Tii  J*([ ^/Otff]b)  *  Tit 


(68) 


na  -  n J0  IWO^r)*  na9  N0 (ttjtkfffr) 


(69) 


Ta  •  T'J0  (l4a/09Hf'r)  *  T\ [fa/Otffr)*  Ta1*  T„  tnr  (to) 

In  these  equations  the  boundary  conditions  at  r  =  0  are  taken  into 
account  but  not  those  at  r  =  R. 


_  dni  dn 
4  a  *  dr  dr 


.£2s  .  r.  -  r„  ,  HL.dTa. 


dr  dr  <71> 


These  conditions  secure  finite  values  of  the  particle  current  densities 
and  finite  values  of  the  sources  at  r  =  £  . 

The  application  of  the  conditions  (71)  at  r  =  f  and  the  conditions 
(65)  at  r  =  R  to  the  equations  (67)  -  (70)  results  in 


".T.fMwf?) 


(72) 


nat  ([*<>/ Oiff]  f)+ ngg  N0  ([*a/ J  *  f) 


16 


(j^  )*  { "a.  J»  (fa/%]*? )♦"„  N,(f“./Ojlf )) 
rh  Jc( !*?)*■  h  • 

f,J.  (fa/%JJf)  *  T"N0  ([<*./%]*?)  i  Ta,  *  T„  tnr 

(%fn,3,(k/Wlt)m 

(£/  { T%  ^MfaA^fJK 

"a,  Jo  (fa/% JlR)'  n«  N»  ffa/%]’  R)  '  ° 

T-w  “  Tqi  *  Ta*  R 

Clearly  a  general  solution  of  this  system  of  equations  is  impracticable. 
We  therefore  restrict  ourselves  to  the  case  of  small  values  of  ?  using 
the  condition 

(78)  i  » 

This  is  not  really  a  new  restriction  since  it  is  included  in  our  assump¬ 
tion  of  a  strong  radial  dependence  of  the  neutral  gas  temperature  as 
may  be  seen  from  our  equation  (66). 

With  equation  (78)  it  is  possible  to  develop  all  special  functions  occu- 
ring  in  equations  (72)  -  (77)  and  to  neglect  higher  order  terms.  With 
equation  (66)  and  equations  (72)  -  (77)  we  have  seven  simultaneous 


(73) 


(74) 


(75) 


(76) 


(77) 
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equation*  for  the  quantities  n^,  n^,  T  j.  T^,  Tq,  J  and  . 

The  quantise*  n^,  n^.,,  T^j,  T^,  Tq  and  f  define  uniquely  the 
electron  ter  ipc mature  ord  density  distributions  in  the  discharge.  The 
quantities  Tq,  f  and  #6ft  are  also  of  interest  with  respect  to  our 
application  in  the  instability  calculation  of  chapter  (VIII). 


With  respect  to  the  application  in  Chapter  (VIII)  we  concentrate  on  the 
calculation  of  the  latter  quantities. 


We  use  the  equations  (66  ),  (  74  )  and  (  76  )  develop  equation  (  66  ) 
and  introduce  it  into  equation  (  76  ).  This  results  in 


j*  lfa/oaffr*)/N.  ([<./%]**) 


7  X 


kTj 


which  together  with 


e(/i4*^.)0^no  •u, 


4*. 

E* 


(79) 


r  .  JL  ♦  r  (so) 

X  *o  '* 

defines  Tq  and  ihg  .  These  equations  contain  apart  from  Tq  and  iff 
the  quantity  E  which  we  express  by  tia  according  to  equation  (  60  ). 
We  see  from  equation  (  60  )  that  the  quantity  E  varies  very  slowly 
with  sC  g  and  therefore  we  use  the  approximation 

*0  "  */•-/>  9Xp('dp/t )  (81) 

where  Eg  is  the  longitudinal  electric  field  oi  the  corresponding  normal 
glow  column. 


Then  the  two  equations  for  T  and  oCg  read 

r°  x  «.  4 


(82) 


and 


v  -<r-  ..-.tv- 
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J.  ([‘‘J04]i*)/N.(fa/0,ll]iR)  * 

<r  » »  .ft 


(83) 

/TV 

kTfa 

mrnJ  tUi  1 

*Ul 

e.  "nTt) 

Tfl  is  the  electron  temperature  in  the  normal  glow  column. 


To  calculate  Tq  and  the  electric  field  E,  respectively  oCa  »  we 
would  have  to  introduce  experimental  data. 

However  without  detailed  calculation  the  following  general  conclusions 
can  already  be  drawn. 

Remembering  that  the  right  hand  side  of  equation  (83)  can  assume  po¬ 
sitive  values  only  and  that  we  should  include  the  transition  to  the  nor¬ 
mal  glow  theory  the  inequality 

<»4>  0,9  it 

must  hold  independent  of  the  choice  of  the  experimental  parameters. 

At  least  this  is  true  within  the  frame  of  our  approximations.  It  follows 
that  the  longitudinal  electric  field  component  E  differs  only  little 
from  the  Schottky  value  Eg.  With  increasing  particle  number  den¬ 
sity  in  the  center  there  is  a  slight  decrease  in  E. 

With  increasing  discharge  current  the  particle  number  density  in  the 
axis  increases  slightly  more  than  proportional  to  the  current. 

Equation  (82)  shows  that  the  axial  electron  temperature  Tq  increases 
with  the  particle  number  density  in  the  axis  and  with  the  total  discharge 
current. 

In  judgement  of  our  results  let  us 'always  remenhber  that  we  must  ex¬ 
clude  small  discharge  currents  due  to  our  initial  assumptions.  Also 
very  large  discharge  currents  are  excluded  since  in  this  region  the 
glow  model  breaks  down. 
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C  GLOW  COLUMN  WITH  NEGATIVE  IONS. 


In  this  chapter  we  consider  the 
influence  of  negative  ions  on  the 
normal  glow  column.  Subject  of 
this  investigation  is  therefore  the 
model  described  in  chapter  (VII  A) 
extended  by  the  addition  of  a  nega¬ 
tive  ion  component. 

Consequently  we  have  the  following 
system  of  equations  describing  our 
problem: 

The  particle  conservation  laws 


o  I 


Pig.  frt  Schematics!  sketch  of  the  density  and 
temperature  distribution  in  the  glow- 
column  with  negative  ions. 


(85) 


(86) 


(87) 


(88) 
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the  neutral  particle  balance  under  the  assumption  of  constant 
pressure,  weak  ionization  and  small  energy  input 


! 


i 


nn  ■  const 


and  the  boundary  conditions 


(90) 


— 1  -  o  for  r  •  o  .  n;  -  o  for  r  •  R 
dr  '  ' 


Introducing  equation  (86)  into  equation  (88)  gives  the  radial  electric 
field 


(91) 


£  _  nh 


and  with  that  the  radial  particle  current  density 


(92) 


r  •  -  D;  Vf  n, 

,r  1  r  *  Zl0kl /iHnH 


In  the  above  equations  the  indices  can  assume  the  values 
(e),  (+)  and  (-)  referring  to  electrons,  positive  and  negative  ions 
respectively.  We  stress  particularly  that  in  this  chapter  -  and  only 
in  this  chapter  -  the  index  (-)  refers  to  the  n  e  g  a  t  i  v  e  ions  and 
not  to  the  electrons. 

For  a  three-component  system  consisting  of  electrons,  singly  charged 
positive  and  negative  ions  the  following  processes  contribute  to  the 
net -production  terms  Aj  ; 


V 

{ 


t 
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Ionization  of  neutral  particles  by  electron  collision 

■  («7p)lp /"„)*«!  n„nt  (93) 

Attachment  of  electrons  to  neutral  particles 

\  -  iiiyp)(p/nx)vtdn'<n*  <«> 

Detachment  of  electrons  from  negative  ions  by  electron  collision 

^  veih  n-ne 

Recombination  of  negative  ions  with  positive  ions 

&d  mCrji  y!th  n>n- 

From  the  equations  (85),  (87),  (92)  and  (93)  -  (96)  we  derive  the  fol¬ 
lowing  simultaneous  differential  equations  for  the  electron  density  rig 
and  the  negative  ion  density  n  . 


(95) 


(96) 


J_  ±  r  f  r'»(j*eO*AO*)+  2m  De  ffeg  7  . 
r  afr  L  n0  +/*)  ♦  2/in.  dr 

n,  [if- -  n-  ] 


(97) 


i  si 

r  dr 


r 


f  Ddn.  (0**0)s*ri-  _ 
i  u  dr  ne(/40 */*)*&". 


**  ** -  T  z*  ytH, 


(98) 
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with  the  boundary  conditions  (90). 


EVALUATION. 


We  do  not  know  of  an  electro  -  ne¬ 
gative  gas  which  can  be  described 
by  a  three -component  system  and 
whose  experimental  data  are  suf¬ 
ficiently  known.  Therefore  we  have 
to  apply  our  calculations  to  a  fic¬ 
ticious  gas.  In  this  model  gas  we  replace  the  various  ion  components 
by  one  positive  and  one  negative  ion  component  with  average  parameters. 
We  will  also  choose  the  numerical  values  close  to  those  known  for  oxygen. 

For  the  evaluation  we  introduce  the  following  abbreviations 

*  mT/*  ,  Ym  ne/nco  ,  z  “  n- 1 n to 

(99) 

*  ‘ l(p*f*«,]/[p°fclT)} 

and  find 


Fig-  7:  Numerical  \  .lines  of  the  experimental  parameters. 
The  right  ordinate  gives  lo'**  lt,/l  resp.  lo* 
v,.<j  '>■"/•> 


(100) 


- till _  *1  S 

*  cfx  y  *2(/t//i€)z  dz  y  L  p  P  0  V»d  z  nn  P  J 


and 


L±)(LLd£t _ £ _  d£  ]  . 

(101)  *  d*  L  T§ch f  y  *  *(/•/&)*  dx  J 

P  I  U4  Vmt/nnPJM  *W  *  nn 


la 

P 
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with  the  boundary  conditions 


e/y 

— —  m  O 

dz 

— —  •  o 

for  x  •  o 

dx  ' 

dx 

/  ’  *  , 

z  •  o 

for  x  •  1 

In  Fig.  7  the  parameters  entering  equations  (100,101)  are  shown  as  a 
function  of  E/p 

The  numerical  values  used  for  the  coefficients  are 

S  «  ro*[cm9]  m/m»  *  to'2 

(103) 

CT  •  1o~n[cm]  T+  •  T  •  Tn  -  Loo  [mKj 

The  system  of  differential  equations  (100,101)  was  solved  by  an  appro¬ 
priate  approximation  procedure  with  the  help  of  a  digital  computer. 

The  results  are  shown  in  FigB.  8-11. 


DISCUSSION. 

The  parameters  of  our  calculation  are  n  /n  aiM*  E/p. 

The  solution  of  the  eigenvalue  problem  produces  the  product  of  the 
total  current  with  the  pressure  ,  Ip  ,  and  the  product  of  the  discharge 
radius  with  the  pressure,  Rp  .  The  corresponding  relations  are  shown 
in  the  Figs.  8  and  9. 
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Fig.  *>i  Kvladvt*  «•!«•*  Iron  «lvii«ilv  n^o  1  nn  •••  ••  iunchon 
of  Iht*  «>MfM*rniHinl<«l  |m ffiiiH'lv m  H|>  •mil  |>l.  (IM  it  thi* 
neuir.il  ga*  «U-n»ily. 


Fl||.  The  reduced  electric  fluid  C/p  ae  a  function  of  the 
experimental  parameter*  kp  md  pi. 


Together  with  the  eigenvalues  Ip,  Rp  the  evaluation  of  the  problem 
produces  a  2 -manifoldness  of  radial  electron  density  distributions. 

We  normalize  these  distributions  using  y  =  na/nan  and  x  =  r/R. 

We  further  introduce  the  relative 
halfwidth  h/R  as  a  characteris¬ 
tic  parameter.  Then  we  find  that 
curves  with  the  same  constriction 
parameter  h/R  show  practically 
the  same  radial  dependence.  Con¬ 
sequently  our  2 -manifoldness  is 
reduced  to  the  1-manifoldness 
shown  in  Fig.  lo.  ri».  101  .Ham,  domlly  dUtrllnitiMu. 

it  tht  constriction  parameter. 

The  electron  density  distributions  may  therefore  be  characterized  by 
the  constriction  parameter  h/R  only.  In  Fig.  lo  we  have  given  this 
typical  parameter  as  a  function  of  the  experimental  data  Rp,  Ip. 

The  range  of  the  experimental  parameters  Rp,  Ip  is  limited  by  the 
assumptions  underlying  our  model.  Since  in  a  collision  dominated  plasma 
the  mean  free  path  of  the  carriers  has  to  be  much  smaller  than  the 


y  =  n  /n  and  x  =  r/R. 
'  e'  eo  ' 
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dimensions  of  the  vessel  it  fol¬ 
lows  that  we  have  a  lower  limit 
for  the  product  Rp,  which  is 
indicated  in  the  dotted  line  in 
the  Figs.  8,9,11.  An  upper 
limit  for  the  product  Ip  is 
required  by  our  assumption  of 
constant  temperature.  This 
upper  limit  is  found  from  a 
careful  consideration  of  the  energy  balance  and  is  approximately 
identical  with  the  largest  value  shown  in  Figs.  8,  9, 11. 

From  the  Figs.  8-11  the  following  conclusions  may  be  drawn: 

A  negative  ion  component  can  cause  remarkable  constriction  of  the 
positive  column. 

For  a  constant  value  of  Rp  the  constriction  increases  with  the  cur¬ 
rent  as  shown  in  Fig.  11.  The  constriction  is  more  pronounced  in  the 
range  of  small  Rp  values  and  practically  negligible  for  large  Rp 
values.  For  constant  current  the  constriction  increases  with  pressure. 

The  characteristic  (E-l)  of  the  column  can  be  found  from  horizontal 
intersections  in  Fig.  9.  It  shows  a  negative  slope.  This  effect  in¬ 
creases  with  decreasing  discharge  radius. 

We  stress  particularly  the  fact  that  the  constriction  effect  calculated 
above  is  nonthermal.  By  this  we  mean  that  within  the  range  of  our 
experimental  parameters  the  assumption  of  constant  particle  compo¬ 
nent  temperatures  across  the  whole  discharge  is  valid.  If  we  would 
proceed  to  larger  values  of  Ip  the  neutral  gas  temperature  develops 
remarkable  radial  dependence  and  with  that  also  the  electron  tempera- 


Fig.  II:  Constriction  parameter  h/R  aa  a  function  of  Use 
experimental  parameters  Rp  and  pi. 
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ture  becomes  a  function  of  the  radial  coordinate.  By  the  same  token 
the  production  coefficient  increases  in  the  center  and  decreases  to¬ 
wards  the  wall.  These  effects  may  cause  an  additional  constriction. 
We  described  the  influence  of  the  thermal  effects  in  the  preceding 
chapter  (VII  B). 

It  is  interesting  to  note  that  already  Gilntherschulze  ,  Seeliger 

34  79 

and  S  o  m  m  e  r  m  e  y  e  r  ’  expected  as  one  cause  of  the  constriction 

34 

process  the  presence  of  negative  ions  and  Gtin  t  h  e  rschulze 
formulated  in  this  connection  the  concept  of  the  "ion  shield  effect". 


D  THE  WALL  DETACHED  GLOW  COLUMN. 


The  preceding  calculations  of 
the  glow  column  are  distinguished 
by  the  underlying  physical  model. 

They  have  in  common  the  asump- 
tion  that  at  the  wall  the  Schottky 
boundary  condition  is  effective. 

We  discussed  already  in  chapter 
(VI)  that  this  boundary  condition 
may  be  rendered  ineffective.  For 
instance  there  can  be  a  competi¬ 
tive  selecting  mechanism  -  like 
for  instance  gas  convection  -  which 
renders  the  simple  particle  balance 

(38)  invalid.  Or  even  the  statistical  fluctuations  of  the  system  may  be 
sufficient  to  outrule  S  chottkyh  condition  if  the  particle  current  den¬ 
sity  -  that  means  the  gradient  of  the  particle  number  density  -  has 
become  sufficiently  small  in  the  neighbourhood  of  the  wall.  In  the  first 
case  the  solution  is  selected  by  the  new  dominating  mechanism  which 
may  be  taken  into  account  by  a  new  effective  boundary  condition.  In 


Fig.  12:  Schematical  skatch  of  ths  density  rad  the  temperature 
distribution  in  the  wall  detached  glow  column.  M|  mo¬ 
dal  sons  of  tbs  normal  glow  approximation.  Mg  modal 
sooa  of  tbs  glow  approximation  with  volume  recombi¬ 
nation  affects. 


the  second  case  the  solution  must  be  determined  from  stability  consi¬ 
derations. 

If  the  boundary  condition  (38)  is  not  effective  any  more  but  the  glow 
approximation  is  still  applicable  then  we  call  the  system  a  wall  de¬ 
tached  glow  column. 

Let  us  consider  a  discharge  volume  in  which  the  glow  approximation 
is  applicable  and  let  us  increase  the  diameter  of  the  discharge  vessel. 
Then  -  keeping  the  total  current  constant  -  the  particle  number  density 
and  the  importance  of  the  divergence  term  will  decrease  and  we  will 
reach  a  limit  where  volume  recombination  and  electron  attachment  be¬ 
come  important.  This  effect  will  first  come  into  play  in  the  outer  re¬ 
gions  of  the  discharge.  It  is  not  our  intention  to  present  here  a  calcu¬ 
lation  of  these  volume  effects.  Fortunately  we  can  refer  the  reader 
to  the  preceding  chapter  describing  the  glow  column  with  negative  ions. 
Here  particle  attachment  and  volume  recombination  is  an  important 
phenomenon.  The  numerical  evaluation  of  the  equations  of  this  chapter 
showed  that  with  the  increasing  contribution  of  these  volume  effects 
the  gradient  of  the  particle  density  decreases  rapidly  near  the  wall. 
This  is  just  the  effect  which  renders  the  Schottky  condition  invalid 
and  causes  wall  detachment. 

E  WALL  STABILIZED  ARC  COLUMN. 

In  the  preceding  investigations  we  applied  the  "glow  approximation". 

The  case  of  the  wall  stabilized  arc  column  is  characterized  by  the  fact 
that  the  "arc  approximation"  described  in  chapter  (VI)  is  applicable 
throughout  the  whole  cylindrical  discharge  volume. 
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Equation  (32)  provides  three  relations  for  the  energy  balance  of  elec¬ 
trons,  ions  and  neutrals. 

Assuming  stationarity,  quasineutrality  and  equal  temperatures 

T  =  T  =  T  the  first  term  on  the  left  hand  side  and  on  the 
+  —  n 

right  hand  side  of  equation  (32)  vanishes.  Multiplication  by  /2 
and  addition  results  in  the  relation 


(104) 


r  dr  (r*  dr  )  f  *  "  "r  dr  (rUifr)  *  5 


Here  the  extensive  quantities  without  index  describe  quantities  summed 
over  all  particle  components.  For  instance  the  quantity  S  represents 
the  radiation  of  the  neutral  particles,  the  electrons  and  the  ions, 
is  the  total  heat  conduction  coefficient  of  all  components.  The  first 
term  on  the  right  hand  side  arises  from  the  ambipolar  diffusion  of 
ionization  energy.  Diffusion  of  excitation  energy  is  neglected. 

We  have  further  from  equation  (31)  the  momentum  balances 

(105)  -  tjx_  n  E  ♦  e  Vn 

(106)  -  •/*♦  n  E  -  e0+  Vn 

The  neutral  particle  momentum  balance  and  conservation  law  will 
be  replaced  by  the  relation 

(107)  nf ,  rn  -  p  •  const 

due  to  the  assumption  of  constant  pressure  and  weak  ionization. 
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We  further  have  the  two  particle  conservation  laws  for  electrons 
and  ions.  However  with  our  assumption  of  quasineutrality  and  the 
basic  concept  of  our  arc  approximation  we  may  use  instead  of  the 
particle  conservation  law  the  charge  conservation  law 


*(},*?-)  *  0 


(108) 


and  the  Saha  equation 


hi 

Further  Maxwell's  equations  provide  the  relation 


(109) 


E  •  "  V  U 


(no) 


The  equation 


7.  £  -  4  7T e  (nf-  />.)  (in) 

defines  the  deviations  from  quasineutrality  and  is  only  of  interest  for 
higher  approximations. 

Altogether  we  have  six  unknown  quantities 

u,  i,J. T  ai2) 

and  the  six  equations  (104-106)  and  (108-110),  assuming  that  the  radiation 
term  S  in  the  energy  balance  is  expressed  by  the  variables  n  and 
T.  The  form  of  this  relation  depends  on  the  state  of  the  plasma,  in 
particular  on  the  question  whether  the  plasma  is  opaque  or  not.  Since 
we  are  considering  a  cylinder  symmetrical  problem  it  follows  from 
(lo8) 
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t- 


(1U; 


where  the  constant  C  must  be  identical  zero  since  the  sum  of  the 
radial  current  densities  should  be  finite  in  the  axis.  By  elimination 
of  the  electric  field  the  momentum  balances  (105,106)  provide  for 
the  radial  particle  current  density  the  relation 


(114) 


n  drt 
}-r  "  e  ^am  dr 


and  for  the  longitudinal  components 


(115) 


•  •/*-  1  '  •/**  n  £* 


Introducing  the  radial  current  density  and  the  longitudinal  current 
density  into  the  equation  (104)  and  remembering  that  the  longitudinal 
component  of  the  electric  field  is  constant  due  to  Maxwell' s 

equations,  we  arrive  at 

(ii6)  r£(r*£)*el/t,*M.)nE',  ' 

~rdr(re  Mi  0am  5?  )  *  £ 

This  equation  (116),  the  Saha  equation  (109)  and  the  function  S  (n,  T) 
expressing  the  radiative  losses,  determine  together  with  the  boundary 
conditions  the  temperature  distribution  in  the  wall  stabilized  thermal 
column. 

As  boundary  conditions  we  have  due  to  symmetrical  reasons  and  accor¬ 
ding  to  equation  (44) 

<117)  dr  m  °  for  r*°  i  Tw  for  r-  R 
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The  equation  (116)  is  still  very  complicated.  A  simplification  can  be 
reached  by  neglecting  the  influence  of  the  ambipolar  diffusion  of  ioni¬ 
zation  energy  in  comparison  to  the  normal  heat  conduction. 

Then  we  have  the  Elenbaas  Heller  differential  equation 

r  Jr  (rX  d?)*e  (>A  7-)n  £*  *  5  (u8) 

The  formulation  of  the  radiative  losses  requires  further  discussion. 

Under  the  assumption  of  a  transparent  channel  we  would  in  principal 
have  to  sum  over  all  emissions  of  all  lines  and  the  continuum.  This 
is  a  too  difficult  problem.  The  usual  procedure  therefore  replaces 
all  levels  by  one  effective  energy  and  writes  the  emission  in  the  form 

(T)  -  (A/T)exp^~  j  (119) 

where  A  is  proportional  to  the  pressure  and  the  average  transi¬ 
tion  probability. 

With  increasing  energy  input  the  contribution  of  the  continuous  radia¬ 
tion  becomes  more  important.  Consequently  the  average  energy  level 
in  our  formula  (119)  rises  . 

If  one  increases  the  energy  input  even  further  one  reaches  the  region 
in  which  the  channel  is  still  transparent  but  the  contribution  of  the 
discrete  spectrum  is  negligible  in  comparison  to  the  contribution  of 
the  continuous  spectrum.  In  this  region  we  can  describe  the  radiative 
loss  by 


S(T)  » 


(120) 


V-- 


with 

n* 

(121)  ^ ^  *  *• 36  '  to****  (kT)i 

and  h  •  'ig  characterizing  the  energy  interval  near  the  ionization 
boundary  through  which  the  energy  levels  may  be  considered  "smeared 
out". 


Depending  on  the  situation  we  introduce  either  (12o)  or  (119)  into  (118). 
Then  we  have  together  with  the  boundary  conditions  (117)  an  eigenvalue 
problem  which  provides  the  temperature  distribution  and  the  correspon¬ 
ding  eigenvalues  of  the  longitudinal  electric  field  E^. 

One  can  reduce  this  eigenvalue  problem  to  an  initial  value  problem 
by  using  the  coordinate  transformation 


which  leads  to  an  equation  of  the  form 

with  the  boundary  conditions 

<i24>  (c<yf  sf  ■  0  f°r  y 0  >  y»  * t 


Here  the  symbols  C^,  C^,  f^  and  f.,  have  the  following  meaning 


c, 


Ak 

AfeUi)* 


Now,  assuming  an  initial  value  of  'v®  and  choosing  the  ratio 
C 2  /  Cl  we  can  integrate  down  to  the  value  ■  0  and  find  a 

corresponding  value  of  y  .  This  value  defines  the  value  of  the 


constant  through  equation  (124). 

the  ratio  of  C 2  /  Cj  both  constants  < 
these  two  relations  we  calculate 
with  equation  (125)  the  experimen¬ 
tal  parameters  Rp  and  .  The 
current  can  be  calculated  from  the 
temperature  distribution  using  the 
Saha  equation  and  the  mobilities 
for  the  electrons  and  ions. 

The  considerations  sketched  above 

have  been  carried  through  by  Wei  - 

100 

zel  and  Schmitz  and 
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K  o  ch  ,  L  e  s  e  mann  ,  W  a  1 1  h  e  r 
Some  of  the  typical  results  are  shown 


Since  we  have  already  chosen 
j  and  C2  are  defined.  From 


Fill.  IJt  Typical  temperature  distributions  fur 
th*  wall  stabilised  arc  [77]  .  The  mea¬ 
ning  of  the  symbols  is  enplainod  in  the 
tout. 


in  Fig.  (13). 


Moreover  it  is  possible  to  define  from  equation  (123)  similarity  laws 
and  draw  general  conclusions  with  respect  to  the  dependence  of  the 
discharge  qualities  on  the  experimental  parameters.  Within  the  frame 
of  our  investigations  here  we  cannot  consider  these  details. 


One  more  word  about  the  radiative  loss.  If  the  channel  of  the  wall-sta- 
lilized  arc  column  is  not  transparent  then  the  theoretical  description 
of  the  radiation  becomes  extremely  complicated  due  to  the  absorption 
and  reemission  processes.  We  do  not  know  of  any  theoretical  treatment 
of  this  problem. 
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F  WALL  DETACHED  ARC  COLUMN. 

In  all  modes  treated  so  far  either  the  "glow  model"  or  the"arc  model" 
was  applicable  throughout  the  whole  discharge  volume.  Only  a  very 
small  region  near  the  walls  -  order  of  magnitude  of  the  mean  free  path  - 
was  excluded.  The  qualities  of  the  "sheath"  were  only  of  interest  for 
the  boundary  conditions. 

The  wall  detached  arc  column  mode  -which  we  investigate  here  -  is 
characterized  by  the  fact  that  both  -  the  arc  and  the  glow  model  - 
are  of  importance. 

To  understand  under  what  conditions  this  wall  detached  arc  mode  occurs 
let  us  consider  the  following  picture. 

We  chose  a  value  of  the  axial  temperature  Tq  and  the  eigenvalue 

E  .  Since  per  definition  the  arc  model  holds  near  the  axis  we  can 
z 

find  the  temperature  distribution 
by  integration  of  equation  (118) 
starting  with  the  boundary  value 
dT/dr  =  o  and  T  =  Tq  at  r  =  o. 

The  corresponding  density  is 
given  by  the  Saha  equation  (109). 

This  integration  may  be  carried 

out  down  to  a  limiting  value  r  =  rc 

at  which  our  arc  approximation 

breaks  down.  The  value  of  r  can 

c 

be  found  from  the  condition  (  132  ). 

We  now  discuss  the  possible  solu¬ 
tions  which  can  exist  for  the  summed 
eigenvalue  and  the  axial  tempe¬ 

rature  T  . 

o 


rig.  Schematic*!  sketch  el  the  deadly  and 
the  temperature  distribution  in  the 
wall  detached  arc  celumn,  rc  tape* 
ratee  the  model  region  ef  the  arc  ap- 
p  real  mat  ion  (Ma)  from  that  of  the  glow 
approximation  (Mg). 
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If  R  *  then  the  arc  approximation  holds  throughout  the  whole 

discharge  vessel.  We  have  a  wall  stabilized  arc  column.  A  solution 
exists  if  the  temperature  of  the  wall  is  identical  with  the  tempe¬ 

rature  calculated  at  the  distance  r  =  R.  T  =  T(R). 

w 

If  the  discharge  radius  R  is  larger  than  r^  our  arc  model  breaks 

down  at  r  and  a  different  model  must  be  used  in  the  outer  part  of 
c 

the  plasma  column. 

This  is  the  glow  approximation  which  is  complimentary  to  the  arc 
approximation. 

Naturally  in  the  environment  of  r  =  r^  there  is  a  transition  region 
in  which  none  of  the  two  models  is  applicable.  The  general  difficulties 
of  such  transition  regions  have  been  discussed  in  detail  in  FTR  1.  In 
accordance  with  the  general  use  we  replace  here  also  the  transition 
region  by  appropriate  boundary  conditions.  To  secure  finite  current 
densities  and  finite  sources  we  have  the  necessary  conditions 


n 


a 


'9  / 


dng 

dr 


±3  T  -  T„  -  SUM 

dr  J  a  9  >  dr  dr 


(127) 


at  r  =  r  . 

c 

The  connection  of  the  arc  approximation  with  the  glow  approximation 

via  equation  (127)  defines  together  with  the  boundary  condition  (39)  one 

and  only  one  value  for  the  discharge  diameter  2R  in  the  range  R  >  r  . 

c 

Our  conclusion  is  therefore: 


If  R  <  r  holds  we  have  a  solution  with  T  and  E  for  any  value 
c  o  z 

of  R  assuming  that  we  can  still  choose  T^.  Should  however  T^  be 
prescribed  then  there  is  only  one  or  no  possible  solution. 
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If  R  >  x  then  the  boundary  condition  (39)  allows  only  one  value  of 
c 

R  for  the  given  values  Tq  and  E^.  However  this  latter  statement 
should  be  carefully  used  remembering  the  discussion  of  chapter  (VI  c). 
There  it  was  stated  that  for  small  gradients  dn/dr  condition  (39) 
may  be  rendered  ineffective.  Other  influences  may  become  effective 
and  determine  R  or  -  if  statistical  fluctuations  prevail  -  all  values 
of  R  are  possible  and  stability  criteria  select  the  observed  mode. 

Under  these  aspects  it  is  of  particular  interest  to  see  if  dn/dr  de¬ 
creases  for  the  wall  detached  arc  column  and  if  the  solution  for  the 
thermal  nucleus  is  only  little  influenced  by  the  position  of  the  wall. 

The  preceding  general  considerations  are  formulated  in  the  following 
calculation: 


For  the  calculation  of  the  thermal  nucleus  we  use  a  simplified  form 
of  equation  (118).  We  neglect  the  radiative  losses.  The  temperature 
dependence  of  the  production  term  via  the  Saha  equation  suggests 
-  as  in  the  case  of  the  normal  glow  column  -  the  representation  by 
a  step  function.  The  position  of  the  step  is  characterized  by  r^.  In 
the  range  o  <  r  <  r^  the  production  term  is  constant.  In  the  range 
r^  <  r  <  r^  the  production  term  is  zero. 

These  are  serious  simplifications  but  they  do  not  affect  the  general 
conclusions  we  want  to  draw. 


The  quantity 
(128) 


r^  is  defined  by 


AT  m  #£*(?,  $ir.)  K  r.* 
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For  reasons  of  finite  energy  current  density  and  finite  energy  sources 
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the  steadiness  of  the  temperature  and  the  temperature  gradient  dT/dr 
must  be  required.  The  the  temperature  diatribution 

T -  7 rc(f-*(%)*}  hr 

r-  r4  {f  -ci(m  A  CJJ  for  r£.r,  <12’ 

follows . 


The  Saha  equation  determines  the  particle  density 

»  -  forrar, 

n  •  K  !('•**•%)]}  f°r  r  2(j 


From  our  criterion  for  the  validity  of  the  arc  approach  (  45  )  the 
critical  density 


(131) 


defines  the  radius  r  to  be 
c 


'c 


r<  arp{h[Tw/K(Kfn’)'t'sJi 


(132) 


In  the  region  r  7  r£  where  the  glow  approximation  applies  we  remem* 
ber  that  the  specific  particle  production  is  small  in  comparison  to  the 
net  production  in  the  arc  center  of  our  discharge.  Consequently  we  neg¬ 
lect  this  contribution.  Then  we  have,  using  condition  (127)  in  the  region 

r  >  r  ,  the  solution 
c 
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(133)  n  •  ne  ( 1  -  J  *•  £  ] 

Condition  (39)  determines  the  value  of  the  discharge  radius  R  . 

(134)  /?■  r,  S)]J 

Two  quantities  are  of  particular  interest.  The  one  is  the  derivative 
dR/dT^  given  by 

(1 33)  *  *  fj  Tu)  <.*'(><*  n»)  *  •} 

and  the  other  the  derivative  (dn/dr)  given  by 

XV 

dn  I  nc  £  ell; 

(136>  dr  I  R  kT0 

r»R 

From  equations  (135)  and  (136)  certain  interesting  conclusions  can  be 
drawn. 


Let  us  first  assume  that  there  are  no  effects  like  convection,  volume 
recombination  or  external  fields  which  could  render  the  Schottky 
condition  invalid.  In  this  case  the  wall  is  the  stabilizing  influence  and 
the  discharge  is  definitely  symmetrical  to  the  axis  of  our  discharge 
vessel.  Equation  (135)  shows  that  under  these  conditions  the  axial  tem¬ 
perature  Tq  and  with  that  the  profile  of  the  thermal  nucleus  is  only 
weak  dependent  on  the  diameter  2R  of  our  discharge  vessel.  Since 
the  thermal  nucleus  determines  the  observed  phenomenon,  this  pheno¬ 
menon  is  practically  independent  of  the  size  of  the  discharge  vessel. 

Then  from  equation  (136)  we  can  see  that  the  gradient  (dn/dr) 

XV 

is  small.  Since  in  many  discharges  we  have  additional  influences  like 
volume  recombination,  convection  etc.  ,  we  can  expect  that  these  effects 
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render  the  S  chottky  condition  (39)  invalid.  Therefore  the  thermal 
nucleus  of  the  discharge  is  not  forced  any  more  to  occupy  the  center 
of  our  discharge  vessel  but  under  the  influence  of  other  effects  may 
show  a  different  geometrical  shape. 

Even  if  the  latter  case  occurs  the  conclusion  that  the  axial  tempera¬ 
ture  and  the  temperature  profile  is  little  influenced  by  the  parameter 
R  will  still  hold  in  a  good  approximation. 

In  summarizing: 

Our  calculations  show  the  existence  of  the  "wall  detached  arc  mode". 

If  the  boundary  condition  (39)  at  the  wall  governs  the  solution  we 
find  that  the  discharge  body  is  axially  symmetric  and  that  the  maxi¬ 
mum  temperature  Tq  is  rather  independent  of  the  diameter  of  the 
discharge  vessel  2R. 

We  further  recognize  that  for  the  wall  detached  arc  mode  with  increa¬ 
sing  R  the  boundary  condition  (39)  may  be  outruled  by  other  influ¬ 
ences  like  convection,  external  fields ,  volume  recombination  etc.  In 
this  latter  case  the  geometry  of  the  wall  detached  column  is  not  defined 
anymore  by  the  presence  of  the  wall  but  by  the  qualities  of  the  stabili¬ 
zing  influence.  If  the  stabilizing  influence  shows  assymmetric  features 
then  the  geometry  of  the  wall  detached  arc  column  will  be  assymmetric. 
In  this  latter  case  too  the  axial  temperature  and  the  thermal  profile 
of  the  nucleus  of  the  wall  detached  arc  column  should  be  practically 
independent  of  the  effective  diameter  of  the  stabilizing  influence. 

It  is  often  quite  difficult  to  relate  experimental  observations  to  theore¬ 
tical  modes.  Experimental  and  theoretical  knowledge  may  be  incom¬ 
plete.  Also  there  are  frequently  non-typical  cases. 
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For  the  wall  detached  arc  column  however  there  are  obviously  cor¬ 
responding  observations. 

78 

Seeliger  described  for  the  first  time  the  continous  change 
from  the  wall  stabilized  glow  to  the  wall  detached  glow  and  the  dis- 
continous  change  to  a  mode  which  he  called  "Schlauchentladung"  and 
which  might  correspond  to  our  wall  detached  arc  column.  Similar 
observations  with  careful  experimental  study  have  been  carried  out 
byEdels,  Gambling  and  Price^’^'^’^’  .  They 

find  that  the  highly  contracted  -  our  wall  detached  arc  column  - 
appears  at  a  certain  characteristic  gas  temperature. 

The  most  recent  results  follow  irom  the  very  interesting  investigation 
of  Kenty  45,46,47  ££n(js  tjiat  ander  certain  conditions  his 

diffuse  "normal  Schottky  mode"  changes  discontinously  to  a  very 
narrow  "filamentary  discharge".  This  filamentary  discharge  shows 
all  the  qualities  which  we  have  described  here  for  the  wall  detached 
arc  column. 

(A  temperature  difference  for  electrons  and  ions  is  not  principally 
excluded  in  our  arc  approximation  or  for  our  wall  detached  arc  column). 


rt«.  ihlMaiA  dtsefcirf  Hi  1mA  «lhN  «Mlsr|s  M  M 


X*  (41). 
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The  investigation*  of  Kenty  prove  that  the  thermal  nucleus  of  the 
wall  detached  discharge  may  be  displaced  within  the  container  without 
difficulties  by  very  small  forces.  The  position,  the  length  and  the  form 
of  the  discharge  is  a  very  uncertain  one.  Kenty  has  shown  that  convec¬ 
tive  forces  are  able  to  elongate  the  discharge  column  extensively  and 
even  to  push  it  into  small  extensions  of  the  discharge  vessel.  Also  he 
observes  that  the  discharge  column  forms  a  little  arch  at  the  side  of  the 
cathode  due  to  the  upward  convection  of  the  gas. 


G  CATHODIC  PART  OF  THE  ARC  DISCHARGE. 

In  the  calculations  of  the  s -modes  of  the  column  presented  in  the  pre¬ 
ceding  chapters  the  governing  boundary  condition  defined  the  symmetry 
of  the  model  regions.  For  the  cylindrical  column  therefore  all  model 
zones  showed  cylinder  symmetry. 

Only  in  the  case  of  the  wall  detached  arc -mode  we  pointed  out  that  if 
the  stabilizing  influence  of  the  wall  is  outruled,  the  geometrical  shape 
of  the  thermal  nucleus  of  the  wall  detached  arc  column  is  not  defined 
by  the  geometry  of  the  wall  and  may  show  various  forms.  Which  of 
the  possible  geometrical  forms  we  will  observe  depends  on  the  presence 
of  another  stabilizing  influence  like  e.  g.  external  fields,  convection 
etc.  or  on  the  stability  qualities  of  the  s-mode.  In  this  respect  we  refer 
to  the  investigations  in  chapter  (VI  C)  and  (X). 

In  more  general  terms  we  can  say  that  in  the  case  of  the  wall  detached 
arc  column  the  discharge  volume  does  not  prescribe  the  interface  se¬ 
parating  the  two  model  regions  of  the  "arc  approximation"  and  the 
"glow  approximation".  We  rather  have  a  whole  continuum  of  stationary 
modes  of  different  geometrical  form.  Which  of  them  is  realised  in  a 
particular  experimental  setup  depends  on  influences  not  included  in  the 
calculations  of  chapter  (VII  F).  Either  we  have  strong  additional  stabi¬ 
lising  factors  or  the  observed  mode  is  selected  by  stability  principles. 
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We  are  confronted  with  exactly  the  same  problem  in  studying  the 
poeeible  ■•mode*  in  front  of  the  electrodes  of  an  arc  discharge. 

Here  the  existence  of  more  than  one  model  region  with  the  uncertain¬ 
ty  of  the  interface  is  not  the  exception  -  as  it  was  in  the  case  of  the 
wall  detached  arc  column  >  but  rather  the  normal  form. 

Approaching  the  electrode,  the  conditions  change  from  those  pre¬ 
sented  in  the  undisturbed  column  into  those  of  the  inertia  limited 
region  immediately  in  front  of  the  electrode.  Consequently  we  have 
in  the  electrode  region  always  several  model  regions  seperated  by 
interfaces. 

The  geometrical  form  of  these  interfaces  are  in  principal  again  defined 
by  the  boundary  conditions.  However  we  have  now  two  competitive  con¬ 
ditions,  one  from  the  electrodes  and  the  other  from  the  walls.  In  cer¬ 
tain  regions  they  counteract  each  other  and  reduce  the  stabilizing  ef¬ 
fect  so  that  statistical  fluctuations  can  render  them  invalid.  This 
leaves  the  geometrical  form  of  the  model  undefined  and  allows  the 
existence  of  a  whole  manifoldness  of  s-modes. 

In  an  earlier  treatment  of  the  electrode  components  of  the  arc  discharge 
FTR  I1  we  have  distinguished  three  model  regions  in  front  of  the  cathode. 

The  first  region  was  the  "contraction  region"  which  extends  from  the 
undisturbed  column  down  to  the  beginning  of  the  inertia  limited  region 
in  front  of  the  cathode.  Within  the  contraction  region  the  arc  approxi¬ 
mation  holds.  It  is  similar  to  the  nucleus  of  the  wall  detached  arc  column. 
As  the  name  suggests  the  contraction  region  is  not  of  cylindrical  shape 
but  describes  the  area  in  which  the  diameter  (2R)  of  the  column  con¬ 
stricts  to  the  value  (2Re)  in  front  of  the  electrode.  The  geometrical 
shape  of  this  contraction  region  as  well  as  the  parameter  of  the  end 

contraction  R  are  not  prescribed  by  the  boundary  conditions  due 
© 

to  the  arguments  stated  above. 
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The  second  model  region  is  the  inertia  limited  zone  in  front  of  the 
cathode  in  which  the  charged  carriers  move  according  to  the  laws 
of  the  free  fall. 

Finally  the  third  model  region  is  the  whole  rest  of  the  discharge 
volume  between  the  interface  of  the  contraction  and  inertia  limited 
region  on  the  one  hand  and  the  walls  and  the  electrode  surface  on  the 
other  hand. 

In  the  preceding  report  ^  we  have  discussed  in  detail  the  qualities 
of  the  contraction  region  and  the  inertia  limited  zone.  This  can  be  done 
without  specifying  the  laws  within  the  third  regions.  The  details  of 
the  calculations  are  not  of  interest  here.  But  there  is  one  essential 
result  which  relates  to  our  present  consideration. 

We  stated  above  that  we  expect  in  front  of  the  electrode  a  whole  mani¬ 
foldness  of  s-modes  corresponding  to  the  possible  geometrical  shape 
of  the  interface  between  the  contraction  region  and  the  external  region. 

Since  the  conditions  immediately  in  front  of  the  cathode  are  not  deter¬ 
mined  and  influenced  by  the  cylindrical  walls  a  circular  onset  might 
in  principal  be  just  as  possible  as  a  quadratic  or  an  elliptic  onset. 

Moreover  we  could  imagine  that  a  stationary  mode  exists  with  more 
than  one  circular  cathode  spot  or  even  with  cathode  spots  of  different 
values  of  end  contraction  meaning  that  the  cathode  spot  area  might  be 
characterized  by  a  set  of  end  contraction  parameters  l/R«t  • 

Furthermore  the  shape  of  the  interface  connecting  the  column  with  the 
cathodic  onset  could  show  various  geometrical  forms. 

i 

Our  investigation  in  FTR  I  pointed  out  that  this  latter  influence  of  the 
geometrical  form  of  the  contraction  region  is  of  little  importance  for 
the  stationarity  condition  in  the  cathodic  and  anodic  area  of  an  arc  dis¬ 
charge.  i 

i 

\ 
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The  situation  is  however  quite  different  with  respect  to  the  pattern 
of  the  end  contraction.  Unfortunately  the  calculations  are  so  extremely 
involved  that  we  could  follow  up  only  the  case  of  a  spherical  single  type 
cathode  onset.  Therefore  we  have  only  one  parameter  of  end  contrac¬ 
tion  l/Rg.  We  investigated  the  possibility  of  stationary  modes  within 
the  range  1  ^  R/R^.  *  .  We  found  that  depending  on  the  experimen¬ 

tal  parameters  stationary  modes  may  be  expected  in  all  ranges  of  the 
parameter  l/R^  and  we  distinguished  three  different  modes  according 
to  the  physical  mechanism  governing  the  current  continuity  in  the 
i.nertia  limited  region. 

There  is  a  thermal  spot  mode.  In  general  the  thermal  mode  occurs  in 
the  range  R/R^  a  *  •  Under  these  conditions  the  energy  transport 

from  the  arc  to  the  electrode  is  large.  It  is  based  on  heat  conduction, 
ion  bombardment  and  radiation.  The  electrode  is  on  high  temperature. 
The  gas  temperature  differs  only  little  from  the  electrode  temperature. 
The  process  governing  the  condition  of  current  continuity  in  the  inertia 
limited  region  is  the  thermal  electron  emission  from  the  cathode  surface. 
Due  to  the  high  temperature  and  the  moderate  current  densities  this 
process  is  able  to  provide  the  necessary  particle  current  component. 

For  the  same  reason  the  influence  of  field  emission  or  more  exactly 
temperature  field  emission  (T-F  emission  FIRI)  is  only  of  little  im¬ 
portance.  Also  the  ion  component,  flow  from  the  gas  to  the  cathode 
surface  is  small  due  to  the  small  gas  temperature. 

The  second  mode  is  the  "combined  spot  mode  (II)".  This  mode  occurs 
in  the  range  of  medium  values  of  end  contraction.  It  is  characterized 
by  the  fact  that  the  combined  action  of  electrons  emitted  from  the  ca¬ 
thode  as  well  as  ions  coming  from  the  gas  governs  the  current  trans¬ 
port  in  the  inertia  limited  region.  The  electrons  may  be  deliberated 
from  the  cathode  by  the  temperature  effect  or  by  the  influence  of  the 
electric  field.  Also  it  is  quite  possible  that  the  process  of  individual 
field  emission  is  of  importance.  Due  to  the  increase  of  the  contrac- 
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tion  parameter  R/R  the  gas  temperature  in  front  of  the  electrode 
can  be  much  higher  than  the  temperature  at  the  electrode  eurface  and 
results  in  a  markable  ion  component  flowing  from  the  gas  to  the  cathode 
surface. 

The  third  mode  is  the  "field  mode  (III)".  This  mode  occurs  with  extreme 
end  contractions.  The  electrode  temperature  is  small.  Due  to  the  high 
current  densities  a  strong  electric  field  at  the  cathode  surface  exists. 
The  process  of  field  emission  -  in  particular  the  I-F  emission  -  is 
operating  and  provides  large  electron  emission  currents.  The  gas  tem¬ 
perature  in  front  of  the  electrode  may  be  high  but  the  ion  current  com¬ 
ponent  is  nevertheless  small  since  the  contact  area  is  so  extremely 
small. 


It  is  necessary  to  distinguish  two  different  types  (Ilia)  and  (Illb)  of  the 
"field  mode". 


In  the  case  Ilia  the  field  at  the  cathode  surface  is  produced  by  ions 
which  originate  at  least  partly  from  thermal  means.  With  other  words 
the  ion  saturation  current  flowing  from  the  hot  gas  region  in  front  of 
the  cathode  to  the  cathode  surface  provides  an  essential  part  of  the  ion 
current  necessary  to  build  up  the  cathodic  field.  Such  a  process  is 
according  to  the  calculations  in  FTR  I  definitely  bound  to  the  condition 
that  the  radius  of  end  contraction  is  much  larger  than  the  mean  free 
path  of  the  ions  A+ 

For  the  mode  Illb  the  thermal  ion  component  flowing  from  the  gas  to 
the  cathode  surface  is  altogether  negligible.  In  this  case  the  ions 
necessary  to  produce  the  electric  field  at  the  cathode  must  be  created 
by  the  electrons  coming  from  the  cathode.  Unfortunately  with  the  small 
cathode  drop  and  the  corresponding  small  ionization  ability  of  the  elec¬ 
trons  the  number  of  ions  produced  by  a  single  electron  is  very  small. 


In  fact  we  require  -  due  to  the  energy  balance  -  multiple  ionization 

processes  to  produce  an  ion.  Consequently  the  ratio  of  the  electron 

current  to  the  ion  current  density  j_  /  j  must  be  a  very  large  one 

if  the  discharge  is  to  operate  in  the  III b  mode.  This  field  mode  Illb 

is  what  was  originally  understood  by  the  field  emission  process  at 

the  cathode  of  an  arc  and  was  already  criticized  by  Compton. 

One  can  see  from  Fig.  2  of  FTR  I  that  it  requires  current  densities 

9  2 

of  extreme  values  (order  of  magnitude  10  amp/cm  )  or  the  presence 
of  processes  favouring  electron  emission. 

We  emphasize  that  the  various  modes  are  possible  depending  on  the 
experimental  parameter.  The  essential  result  of  our  calculation  in 
FTR  I  was  however  that  for  a  fixed  set  of  experimental  parameters 
the  number  of  possible  stationary  modes  -  characterized  by  our 
contraction  parameter  R/R^  *  was  quite  restricted.  The  solution  of 
the  stationary  equations  of  our  problem  did  not  allow  all  possible 
values  of  end  contraction  R/R^,  but  only  a  certain  set  which  could  be 
determined  by  the  so-called  "EXISTENCE  DIAGRAM".  Such  a  diagram 
defines  the  possible  stationary  solutions  from  the  intersection  of  two 
curves  representing  the  ion  saturation  current  and  the  ion  defect  cur¬ 
rent.  The  points  of  intersection  are  called  POINTS  OF  EXISTENCE 
(E-POINTS)  and  determine  those  values  of  the  parameter  R/Re  which 
correspond  to  the  possible  stationary  modes.  In  general  such  an  existence 
diagram  can  have  up  to  four  intersection  points.  But  there  may  be  also 
one  or  even  no  intersection  points  within  the  range  of  the  stationary 
modes  I,  II,  and  III. 

This  result  means  that  we  do  not  have  -  as  one  might  expect  -  a  con- 

tinous  manifoldness  of  s-modes  for  the  cathode  spot  onset.  The  result 

also  explains  the  phenomenon  of  constriction  of  the  arc  in  front  of  the 

cathode.  It  explains  moreover  the  existence  of  certain  critical  pressure 

and  current  values  (p  ,  I  )  in  very  satisfactory  agreement  with  ex- 

c  c 

perimental  observations. 
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In  conclusion  then  we  state  that  in  front  of  the  cathode  of  an  arc  dis¬ 
charge  we  have  in  general  a  finite  number  of  possible  s-modes.  For 
the  case  of  a  circular  single  type  onset  we  have  four  or  less  than  four 
of  such  modes.  Whether  other  multiple  spot  modes  or  modes  of  diffe¬ 
rent  geometrical  shape  may  exist  in  front  of  the  cathode  hitherto  has 
not  been  investigated.  Buttheie  is  no  doubt  that  other  modes  of  that 
type  exist.  To  prove  this  we  show  in  Fig.  16  pictures  of  several  modes 
of  high  pressure  mercury  arc  discharge  which  exist  for  exactly  the 
same  values  of  external  parameters  and  were  observed  to  change  from 
one  into  the  other  without  any  change  in  the  external  parameters. 


Fi|.  16s  examples  for  ths  different  spot  inodes. 
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H  ANODIC  PART  OF  THE  ARC  DISCHARGE. 


It  is  one  of  the  satisfactory  features  of  our  theoretical  picture  about 
the  electrodic  parts  of  an  arc  discharge  that  it  does  not  underlie 
different  concepts  for  the  cathodic  and  the  anodic  region. 

The  basic  mechanism  which  uses  only  the  well  known  classical  phy¬ 
sical  laws  applies  as  well  to  the  anodic  region  as  to  the  cathodic 
region.  The  difference  observed  in  front  of  the  anode  in  comparison 
to  the  cathode  can  be  caused  only  by  quantitative  differences  in  the 
parameters  entering  the  calculations.  There  are  essentially  two  of 
such  features. 

One  is  the  lack  of  an  anodic  emission  process  corresponding  to  the 
emission  from  the  cathode.  (We  restrict  here  from  the  special  case 
of  the  Beck  -arc  where  we  have  an  ion  emission  component). 

The  second  feature  is  the  difference  in  the  mass  of  the  particle  com¬ 
ponent  providing  the  saturation  current  flowing  from  the  gas  to  the 
electrode.  In  front  of  the  cathode  this  component  is  carried  by  the 
ions.  In  front  of  the  anode  by  the  electrons.  For  corresponding  con¬ 
ditions  in  the  end  contraction  region  the  saturation  currents  flowing 
will  show  the  ratio 

(137)  I./ If’  ("*/'».)* 

It  is  possible  to  explain  the  principal  differences  between  the  cathodic 
and  anodic  parts  of  the  arc  with  the  help  of  this  simple  relation.  We 
refer  to  the  difference  in  the  end  contraction  at  the  two  electrodes  and 
the  so-called  micro-spot  mode. 

Concerning  the  details  of  the  calculations  for  the  anode  onset  we  refer 
to  FTR  I.  The  basic  result  is  again  -  as  in  front  of  the  cathode  -  a 
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limited  number  of  a  •modes.  For  the  single  circular  type  onset  we 
have  four  or  less  than  four  possible  modes.  Whether  for  different  onset 
pattern  of  multiple  structure  and  noncircular  shape  additional  s- modes 
exist  has  never  been  investigated  theoretically  but  as  in  the  case  of  the 
cathode  onset  there  is  experimental  evidence  for  this  fact.  (FTR  I). 


I  CATHODIC  PART  OF  THE  GLOW  DISCHARGE. 

The  "glow  approximation"  and  the  "arc  approximation"  described  in 
chapter  VI  was  used  to  distinguish  the  various  s-modes  of  the  arc  and 
glow  column.  They  are  also  characteristic  for  the  distinction  of  the 
electrodic  parts  of  the  glow  and  arc. 

The  electrodic  parts  of  the  arc  are  described  in  their  main  body  -  the 
contraction  region  -  by  the  assumptions  underlying  the  arc  approxi¬ 
mation.  There  is  a  direct  resemblance  to  the  theoretical  description 
of  the  wall  detached  column.  The  electrodic  parts  of  the  glow  on  the 
other  hand  are  governed  by  the  approximations  characteristic  for  the 
"glow  approach".  At  least  this  is  true  as  far  as  the  particle  conser¬ 
vation  is  concerned.  With  respect  to  the  energy  balance  there  are 
certain  differences. 

We  use  two  experimental  parameters  to  characterize  the  cathodic 
part  of  the  glow  discharge.  One  is  the  extension  1  and  the  other  the 
total  voltage  drop  U^.  We  also  assume  that  we  can  treat  our  problem 
one-dimensionally.  This  latter  assumption  does  not  exclude  a  radial 
limitation  of  our  discharge  area  but  assumes  only  that  the  radial  in¬ 
fluences  are  negligibly  small. 

Our  interest  is  here  focussed  on  the  question  what  stationary  modes 
may  exist  for  a  given  set  of  experimental  parameters  Uc>  1. 
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First  of  all  we  have  to  distinguish  between  two  principly  different 
kinds  of  s -modes  depending  on  whether  the  electric  field  through 
the  distance  1  is  approximately  constant  or  shows  a  distorted 
shape  due  to  the  presence  of  space  charges. 

The  first  case  occurs  in  the  range  of  small  current  densities.  Here 
we  have  only  one  model  region.  The  requirement  of  particle  conser¬ 
vation  has  the  simple  form 


Here  *6  is  the  production  coefficient  per  cm  and  no  volume  recom¬ 
bination  is  assumed. 

The  integration  of  this  equation  gives  the  current  density  of  the 
electrons  at  the  anodic  end  as  a  function  of  the  electron  current 
density  at  the  cathode  surface  in  the  form 

*  * 

(1.39)  i  (Z)  •  }J<>)  «*P{J  ^  J 

"  o 

Due  to  the  current  continuity  the  ion  current  must  be  given  by  the 
formula 

(140)  }+(*)-&<•)  m?J°)[  1m**p{f**(y)d5}] 

0 

If  we  assume  that  no  ions  enter  the  anodic  end  and  use  the  cathodic 
boundary  condition  that  each  ion  at  the  cathode  surface  deliberates  elec¬ 
trons  according  to  the  relation 


(141) 
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then  the  well  known  stationarity  condition  of  the  Towns  end  discharge 
results 

1  '  IT  [**/>{  f*(s)ds  }->]  (142) 

o 

This  is  one  possible  s-mode  for  our  cathodic  region  characterized  by 
U  ,  1.  In  fact  it  does  not  describe  only  one  mode  but  a  whole  series 
of  modes  since  the  parameter  j  does  not  occur  in  the  stationarity 
condition.  Therefore  for  a  discharge  with  a  given  total  current  the 
radial  extension  is  still  a  variable  parameter.  The  condition  is  that 
the  current  density  does  not  exceed  the  limit  prescribed  by  the  space 
charge  effects. 

If  this  last  condition  is  not  met  then  we  approach  a  completely  diffe¬ 
rent  type  of  cathodic  mode.  This  new  type  of  mode  is  characterized 
by  the  fact  that  the  Poisson  equation  causes  a  strong  change  in  the 
field  distribution.  To  describe  this  s-mode  we  must  distinguish  three 
different  model  regions. 

The  first  one  is  the  so-called  cathode  fall  region  (CFR).  In  this  area 
we  find  practically  the  whole  cathode  fall  U^.  The  electrons  move  from 
the  cathode  through  the  cathode  fall  region  producing  in  a  rather  strong 
field  additional  electrons  and  ions.  At  the  end  of  the  CFR  the  electrons 
-  at  least  a  large  portion  of  them  -  is  in  a  "runaway  state  .  We 

stress  that  for  the  existence  of  the  glow  cathode  mechanism  described 
here  the  production  of  electron -runaways  is  decisive.  The  ions  produced 
in  the  CFR  move  towards  the  cathode,  bombard  the  cathode  and  liberate 
electrons.  But  not  only  the  ions  produced  in  the  CFR  reach  the  cathode. 
There  is  also  an  ion  component  entering  the  CFR  from  the  neighbouring 
model  zone.  The  so-called  negative  glow  region  (NGR).  On  the  other 
hand  electrons  coming  from  the  CFR  enter  the  negative  glow. 
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The  second  model  zone  is  the  negative  glow.  Here  the  electrons 
from  the  cathode  fall  have  beam  character  and  enter  the  negative 
glow  region  with  a  certain  initial  energy.  In  passing  through  the 
gas  they  produce  secondary  electrons  and  ions  which  leave  the  nega¬ 
tive  glow  by  ambipolar  diffusion.  Within  the  negative  glow  area  we 
have  hardly  any  electric  field.  Also  the  space  charge  is  negligible. 

As  we  know  from  experiments  there  are  three  energy  groups  in  the 
negative  glow.  The  primary  electrons  which  have  the  properties  of 
a  beam.  The  secondary  electrons  with  an  average  energy  of  order 
or  magnitude  5-10  eV.  These  are  the  particles  produced  by  the 
primaries.  Finally  the  dominating  "ultimate  electrons"  which  have 
low  energy  and  show  practically  a  Maxwellian  velocity  distribu¬ 
tion. 

The  third  model  region  is  the  Faraday  dark  space.  Here  the  produc¬ 
tion  of  particles  is  practically  zero.  The  runaways  from  the  CFR 
have  been  broken  down  in  the  negative  glow.  Only  secondary  particles 
diffuse  through  the  cathodic  end  into  the  Faraday  dark  space.  There 
is  no  carrier  production  and  no  carrier  destruction  in  the  Faraday 
dark  space.  All  particles  diffuse  in  the  axial  and  radial  directions 
and  are  destroyed  by  wall  recombination.  There  is  little  electric 
field. 

To  find  the  stationary  modes  possible  within  the  frame  of  these  model 
regions  we  construct  as  in  the  case  of  the  arc  an  existence  diagram 
plotting  the  extension  1  on  the  abscissa  and  the  voltage  on  the 

ordinate.  (See  Fig.  17). 

The  stationary  states  are  governed  by  two  relations  (FTRII).  One  is  the 
Pois  son  equation  integrated  under  the  assumption  n+  ^  n_  . 

It  gives 
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where  d  is  the  extension  of  the  CFR. 

The  other  is  the  stationarity  law  which  we  find  from  a  calculation 
similar  to  that  for  the  Townsend  discharge  (142)  but  now  limited 
to  the  extension  of  the  CFR.  The  only  difference  to  the  Townsend 
calculation  is  that  we  now  take  into  account  the  possibility  of  an  ion 
component  flowing  from  the  negative  glow  to  the  CFR.  This  provides 
the  boundary  condition 


(d)  •  &  ( d )  (144) 

Under  these  circumstances  we  have  the  stationarity  condition 

1  *  ?  [(***)  m*P  l  J *5 )  '  1  ]  (U5) 

Let  us  draw  attention  to  the  fact  that  our  existence  diagram  here  has 
quite  a  similarity  to  the  existence  diagram  in  arc  theory.  Our  sta¬ 
tionarity  condition  (145)  describes  nothing  else  than  the  current  con¬ 
tinuity  in  our  glow  discharge.  Just  as  the  existence  diagram  in  the 
case  of  the  arc  was  based  on  the  current  continuity  requirement  in  the 
inertia  limited  zone. 

Fig.  17  shows  schematically  A-curves  and  B-curves.  The  A-curves 
correspond  to  the  stationarity  condition  (145).  The  parameter  for  the 
curves  is  the  value  d/l.  The  thin  lines  correspond  to  the  Poisson 
equation  after  eliminating  the  quantity  d  with  the  help  of  the  sta¬ 
tionarity  condition  (145).  The  resulting  law  is  independent  of  1  and 
consequently  is  represented  in  our  -  1  diagram  by  lines  paral¬ 

lel  to  the  abscissa.  The  parameter  of  the  B-curves  is  the  current 
density  j. 
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Any  point  of  intersection 

of  a  B-curve  with  an  A- 

curve  i»  a  poaaible  mode 

(point  of  exiatence)  for 

the  valuea  U  and  1. 
c 

The  parameter  a  d/l,j 
are  defined  by  these  in¬ 
tersection  points.  That 
means  for  a  given  total 
current  J  with  j  the 
radial  exiatence  of  the 
discharge  is  defined. 

Our  existence  diagram  is  quite  revealing.  It  not  only  defines  the  pa¬ 
rameter  values  of  j  and  d.  It  also  shows  that  there  are  areas  in 
which  no  modes  of  the  present  type  exist.  In  the  Fig.  17  we  have  in¬ 
dicated  these  areas  by  I  and  II.  The  existence  of  I  explains  the 
phenomenon  of  a  certain  minimum  current  density  which  is  required 
for  the  cathodic  part  of  the  glow  discharge.  It  also  explains  the  exis¬ 
tence  of  a  certain  minimum  voltage,  which  is  well  known.  The  exis¬ 
tence  of  the  region  II  corresponds  to  the  phenomenon  of  the  obstructed 
glow  discharge.  It  shows  that  below  a  certain  extension  d^  of  the 
cathodic  discharge  parts  the  minimum  voltage  increases  with  decrea¬ 
sing  value  of  1.  This  too  has  been  observed. 

Finally  we  have  indicated  in  our  existence  diagram  Paschen's  law 
which  corresponds  to  the  Townsend  mode  described  above  .  (Aq). 

Our  existence  diagram  gives  the  following  results.  There  is  a  certain 
area  of  the  parameter  values  Uc  -  1  in  which  no  stationary  modes 
of  the  glow  discharge  type  (CFR,  NGR,  FDR)  exist.  Then  there  is  an 
area  of  Uc  -  1  in  which  we  have  always  a  stationary  mode  but  only 
with  a  definite  value  of  the  parameter  d  (extension  of  the  CFR)  and  a 


Vi)!.  ITj  E  •diagram  of  thv  glow  diathargr. 
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definite  value  of  the  current 
density  j  (  radial  extension 
of  the  discharge  for  a  given 
value  of  the  total  current,  J), 

In  the  same  region  there  exists 
for  any  value  of  respective¬ 
ly  1  also  an  additional  mode 
which  has  only  one  model 
region  and  an  undistorted  elec¬ 
tric  field  -  the  T  ownsend 
mode. 


Fill-  IN:  Voltage-  current  density  characteristic  fur  different  values  p^,  in¬ 
cluding  the*  Townsend,  normal  and  abnormal  ulsw  regions. 


So  far  we  have  discussed  the  question  of  stationary  modes  using 

-  1  as  parameters.  It  may  be  desirable  to  consider  the  same 

situation  using  Uc  and  J.  To  find  these  results  we  may  first  take 

from  Fig.  17  the  cathodic  voltage  drop  as  a  function  of  j.  This  is 

demostrated  in  Fig.  18.  As  this 

figure  shows  for  j  w  j  we 

have  an  increase  in  U  with  the 

c 

current  density  which  corresponds 
to  the  anomalous  glow  discharge. 

For  current  densities  j  <■  j 

n 

we  are  in  the  region  of  the 
Townsend  discharge  where  the 
space  charge  is  neglected  and 
consequently  the  voltage  is  inde¬ 
pendent  of  the  current  density. 

Then  however  in  this  region  there 

is  an  influence  of  the  parameter  r‘''  mm i«- 

1  due  to  the  stationarity  condition 
and  consequently  a  splitup  of  the 
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characteristic  in  several  branches.  The  continous  transition  from 
Townsend  to  glow  discharge  shown  in  Fig.  18  is  due  to  space 
charge  effects  and  cannot  be  taken  from  the  calculations  presented 
above.  Again  we  have  the  problem  of  a  transition  region. 

Now  from  Fig.  18  we  may  derive  the  desired  characteristic  U  -  J. 
We  introduce  the  relative  cathodic  area  covered  by  a  discharge  cross 
section  ^  =  F/F^  (Fc  =  cathode  area)  as  a  parameter  and  draw 
the  corresponding  characteristics  as  shown  in  Fig.  19  .  The  para¬ 
meter  1  has  a  lower  limit  determined  mainly  by  the  influence 

of  the  radial  effects,  which  is  indicated  in  Fig.  19  by  the  upward  shift 
of  the  corresponding  (t?'  ~  J)^  characteristics.  Obviously  there  is 
no  unique  relation  between  U  and  J  in  Fig.  19. 

We  may  therefore  conclude: 

There  exists  a  whole  area  Uc  and  J  where  we  ha/e  no  stationary 
solution  of  the  Townsend  or  glow  discharge  type.  Then  there  is 
an  area  in  which  we  have  at  least  two  stationary  modes  of  the  grow 
discharge  type  which  are  distinguished  by  the  parameter  V  des¬ 
cribing  the  coverage  of  the  cathode  surface.  Finally  there  is  a  cer¬ 
tain  number  of  combinations  -  J  where  we  have  not  only  solu¬ 
tions  of  the  glow  discharge  type  but  also  solutions  of  the  Townsend 
type. 
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J  ANODIC  PART  OF  THE  GLOW  DISCHARGE. 

The  situation  of  the  anodic  parts  of  the  glow  discharge  is  quite  diffe¬ 
rent  from  that  of  the  cathodic  parts. 

The  cathodic  part  of  the  glow  discharge  is  of  vital  importance.  The 
anodic  part  is  not.  There  is  a  substantial  theoretical  understanding 
of  the  cathodic  part  and  we  can  describe  various  possible  modes  on 
the  basis  of  one  typical  model,  the  "glow  model".  There  is  hardly 
any  quantitative  theory  for  the  anodic  part  and  there  is  also  no  typi¬ 
cal  anode  model. 

However  the  existence  of  various  s -modes  has  been  observed  though 
they  could  not  be  treated  theoretically.  We  restrict  ourselves  there¬ 
fore  here  to  give  a  qualitative  picture  of  these  observations. 

If  we  first  consider  a  plane  anode  placed  into  the  discharge  and  being 
in  the  negative  glow  then  we  have  no  detectable  light  phenomenon  and 
the  anode  drop  is  negative.  This  effect  has  been  termed  "the  negative 
anode  fall". 

When  the  anode  is  gradually  moved  away  from  the  negative  glow  into 
the  Faraday  dark  space  and  further  into  the  positive  column  then  a 
luminous  layer  forms  in  front  of  the  anode,  -  the  so-called  anode 
glow.  At  the  same  time  the  potential  increases  to  positive  values 
with  respect  to  the  plasma  potential.  These  observations  are  quite 
understandable  on  the  simple  basis  that  the  current  of  the  charge 
carriers  diffusing  to  the  anode  surface  must  be  identical  with  the 
total  current  flowing  in  the  discharge.  Within  the  negative  glow 
the  electron  current  component  is  too  large  and  therefore  a  retarding 
field  is  necessary.  In  the  positive  column  however  the  electron  dif¬ 
fusion  current  is  not  sufficient  and  therefore  an  accelerating  electric 


field  ie  required  which  at  the  same  time  produces  the  light  phenomenon 
of  the  anode  glow.  The  increase  in  the  positive  anode  drop  is  limited 
when  ions  are  produced  within  the  anode  fall  region  so  that  the  positive 
space  charge  can  break  down  the  development.  This  simple  type  of 
anode  mode  can  be  described  roughly  by  an  analysis  which  is  similar 
to  that  applied  to  the  CFR.  We  refer  the  reader  to  FTR  II. 

Modes  different  from  the  simple  picture  given  above  may  be  observed 
in  those  cases  where  the  current  density  in  front  of  the  anode  is  in¬ 
creased  due  to  a  special  geometry  or  limited  surface  of  the  anode. 

Also  if  a  gas  load  is  present  one  may  observe  the  so-called  "pearls 
and  bubbles"  which  present  a  different  anode  s-mode.  Fig.  2o  shows 
some  examples  of  this  type  of  mode.  The  axial  extension  of  these  spots 
is  much  larger  than  the  mean  free  path  of  the  electrons  which  charac¬ 
terize  the  extension  of  the  anode  glow.  The  anode  spots  occur  preferen¬ 
tially  at  corners  and  edges  where  the  inhomogenity  of  the  field  causes 
an  increase  in  the  current  density.  Depending  on  the  experimental  setup 
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there  may  be  a  large  number  of  anode  apota  which  may  be  arranged 
in  a  regular  geometric  pattern.  (Fig.  20a).  If  the  preaaure  ia  de- 
creaaed  the  pearla  develop  into  opaque  bubblea  (Fig.  20b). 

The  model  underlying  the  phenomenon  of  pearla  and  bubbles  is  not  yet 
fully  understood.  Weizel  andHornberg  claim  that  such  a  pearl 
has  an  extension  one  order  of  magnitude  larger  than  that  of  the  normal 
anode  fall  region.  Therefore  the  concepts  of  mobility  and  diffusion 
should  be  appropriate  for  describing  the  situation  in  the  pearls  and 
bubbles.  Weizel  and  H  o  r  nb  e  r  g  used  these  concepts  in  the  P oi  s  -  - 
son  equations  and  at  the  same  time  applied  the  stationarity  condition 
similar  to  the  calculation  for  the  CFR.  In  agreement  with  the  geome¬ 
trical  shape  they  applied  a  spherical  polar  coordinate  system  and  came 
to  the  conclusion  that  the  pearls  may  be  described  as  a  plasma  with 
a  positive  space  charge.  The  plasma  is  an  extensive  one  which  produces 
many  more  carriers  than  are  required  for  the  current  continuity,  into 
the  column.  The  region  of  the  pearls  immediately  in  front  of  the  anode 
is  described  in  linear  form  and  essentially  based  on  the  concept  of 
ambipolar  diffusion  towards  the  anode.  From  these  considerations 
Weizel  and  Hornberg  found  the  potential  distribution.  Again  we 
refer  for  details  to  FTR  II. 

For  our  present  application  the  essential  conclusion  is: 

In  front  of  the  anode  -  as  in  front  of  the  cathode  -  a  limited  number 
of  s -modes  are  possible.  The  normal  planar  form  of  the  anode  fall 
region  can  be  described  in  terms  similar  to  those  of  the  CFR.  In 
addition  "pearl  and  bubble  modes"  have  been  experimentally  observed 
and  according  to  W  e  i  a  e  1  and  Hornberg  may  be  described  as  a 
plasma  with  a  positive  space  charge  and  an  electrodic  sheath  which  is 
governed  by  the  laws  of  ambipolar  diffusion.  The  normal  one-dimen¬ 
sional  AFR-s-mode  and  the  pearl  and  bubble-s -mode  can  exist  in  the 
same  range  of  experimental  parameters. 
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VIII  THERMAL  INSTABILITY  OF  THE  GLOW  COLUMN. 


INTRODUCTION. 


42 

Kadomtsev  and  Nedoipasov  were  the  first  to  tackle  the 

problem  of  the  instability  of  a  collision  dominated  plasma  column  with 
a  superimposed  external  longitudinal  magnetic  field.  In  subsequent 
papers  J  ohn  s  on  and  Jerde,  Holter  and  others  ^ 

treated  the  problem  in  more  detail  and  with  more  mathematical  rigour. 
Sirtion  andKrall  considered  the  Kadomt  s  e  v  type  instabi¬ 

lity  of  a  low  pressure  discharge. 


The  result  of  these  calculations  shows  that  with  increasing  magnetic 

field  the  m  =  1-mode  becomes  unstable  if  the  magnetic  field  surpasses 

a  certain  critical  value  B  . 

c 

The  physical  picture  behind  this  instability  is  rather  simple.  A  density 
disturbance  of  the  stationary  state  produces  an  increase  in  the  particle 
production  on  the  one  hand  and  an  increase  in  the  particle  loss  on  the 
other  hand.  If  the  increase  in  the  particle  production  overcomes  the 
increase  in  the  particle  loss  and  the  two  effects  do  not  change  the  ini¬ 
tial  shape  of  the  disturbance  then  we  have  an  unstable  eigensolution  of 
the  problem. 


One  of  the  basic  assumptions  of  all  the  calculations  quoted  above  is 
that  together  with  the  density  disturbance  no  temperature  disturbance 
occurs. 

Under  these  circumstances  the  mao-mode  is  stable,  independent  of 
the  magnetic  field.  The  onset  of  the  instability  is  determined  by  the 
msl-mode. 
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However  there  are  experiments  which  show  an  instability  of  the 

m=o-mode  with  increasing  current  and  increasing  pressure.  As  an 

47 

example  we  refer  to  the  experiments  of  K  e  n  t  y  and  Hermann 


fig.  Hi  Examples  for  column  modes  before  ami  after  thermal  instability. 


Due  to  the  neglection  of  a  temperature  disturbance  the  calculations  of 
the  Kadomtsev  type  cannot  be  applied  to  these  observations. 

A  temperature  disturbance  may  cause  instability  even  in  the  absence 
of  a  magnetic  field.  It  is  the  aim  of  this  chapter  to  investigate  this 
"THERMAL  INSTABILITY"  of  a  cylindrical  collision  dominated  plasma 
column  with  and  without  a  longitudinal  magnetic  field. 

SYSTEM. 

Subject  of  this  investigation  is  again  an  infinite  cylindrical  plasma 
column  of  radius  R.  The  column  is  homogeneous  in  the  z-direction. 
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It  contains  gas  of  one  kind.  The  theory  of  a  column  with  several 
neutral  particle  components  follows  exactly  the  same  lines  and  re¬ 
quires  only  additional  summation.  The  gas  is  partially  ionized  and 
we  have  electrons  and  one  kind  of  singly  charged  ions.  The  concept 
of  quasineutrality  holds.  A  longitudinal  magnetic  field  of  strength  B 
can  be  superimposed. 

CONCEPT. 

Within  the  range  of  the  normal  glow  column  (Schottky  column)  we 
expect  no  temperature  instability.  As  we  will  see  the  temperature  in¬ 
stability  depends  decisively  on  the  coupling  between  density  disturbance 
and  temperature  disturbance.  In  the  Schottky  concept  the  tempe¬ 
rature  is  assumed  independent  of  the  density  distribution  and  there  is 
consequently  no  coupling  between  density  disturbance  and  temperature 
disturbance. 

The  situation  is  different  in  the  case  of  the  thermal  glow  column  des¬ 
cribed  in  chapter  VII  B.  Here  per  definition  the  energy  balance  is  of 
importance  and  causes  a  coupling  of  the  temperature  distribution  with 
the  density  distribution.  Consequently  a  density  disturbance  by  necessi¬ 
ty  produces  a  temperature  disturbance.  The  temperature  disturbance 
results  in  a  change  of  the  particle  production  coefficient.  In  this  way 
an  additional  influence  on  the  particle  production  term  via  the  electron 
temperature  can  be  expected  and  may  cause  instability  even  if  there 
is  no  instability  of  the  Kadomtsev  type. 

For  instance:  A  spontaneous  increase  in  the  electron  density  near  the 
axis  of  the  discharge  results  in  a  decrease  of  the  neutral  particle  den¬ 
sity  which  in  turn  produces  an  increase  in  the  electron  temperature 
and  with  that  an  additional  particle  production.  This  mechanism  can  be 
the  root  of  an  instability. 


'  -•  -r  v !"-v,~  :■  8/**#  iLIMIIjJU'i'KWQMMMMMMi 
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BASIC  EQUATIONS. 


The  equations  governing  the  thermal  glow  column  can  be  taken  from 
chapter  VII  B.  They  are  represented  in  the  law  of  particle  conser¬ 
vation  of  electrons 


dn  Orf  d_  /_  dn  \ 

JT  r*(rT')-*  n  -  ° 

and  the  equation  for  the  electron  temperature 


(14b) 


“  O 


(147) 


with  the  particle  production  per  second 


*  •  A**P  {-  $} 


(148) 


and  the  coefficients 
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and 
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The  quantities  B  and  AT  are  known  from  experiment  (see  equ.  60). 


Equations  (146)  and  (147)  are  two  simultaneous  differential  equations 
for  the  variables  T  and  n  . 


These  equations  together  with  the  boundary  conditions  in  the  axis 


(151)  ■—  ■  o  ~r~  “  O  for  r  •  o 

dr  •  dr 

and  at  the  wall 

(152)  rt  -  o  f  T-  o  for  r  •  R 

define  the  radial  distributions  of  temperature  and  particle  density 
as  well  as  the  eigenvalue  of  the  longitudinal  electrid  field  E^. 

SOLUTION  OF  THE  EQUATIONS. 


We  want  to  treat  the  stability  problem  by  the  method  of  normal  modes 
which  is  based  on  a  first  order  perturbation  theory  of  the  steady  state 
solution  of  our  system.  We  therefore  use  the  development 

(153)  T-  TM*  Tnl ; 


where  the  index  (o)  refers  to  the  steady  state.  The  index  (1)  designates 
the  disturbance.  Introducing  equation  (153)  into  equations  (146)  and  (147) 
and  remembering  that  the  quantities  with  the  index  zero  satisfy  the 
equations  (55)  and  (63)  we  arrive  at  a  set  of  relations  for  the  first 
order  contributions 


(154) 


and 
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where  all  higher  order  terms  are  neglected. 


We  look  for  the  eigensolutions  of  our  perturbed  problem  by  introducing 
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the  relatione 

n<”  <~,  l)  •  >  f r)  «"■*' ,  r‘V,  fj  •  7>  Cr;*"-*'  (15‘) 

which  gives 
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and 
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Since  the  equations  are  dominated  by  the  operator 


1.  JL  /r  d  \  A 

7  d?<rdF)  m  *r 


(159) 


we  develop  our  eigensolutions  in  a  set  of  Bessel  functions  according 
to 


nM  ,r>  “  |  n*.1  J*  (®»  r/*) 
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where  the  are  the  zeros  of  the  zero  order  Bessel  function  J  . 

o 


From 


| (,'rw  -  ft'/V” 


(161) 


we  find  the  relation 
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(162) 
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If  we  introduce  these  developments  (160)  and  (162)  into  equations  (157) 
and  (158)  a  system  of  simple  uncoupled  linear  equations  for  the  coeffi¬ 
cients  n^(yand  T^  ^would  result  if  the  coefficients  of  the  equations 
(157)  and  (158)  were  all  constants. 


Now  in  our  present  problem  the  coefficients  D,  /*,  »*'  and  X  are 
indeed  considered  constant.  However  the  terms  r/*1  and  ^ 

show  a  strong  radial  dependence  which  is  defined  by  the  solutions  of 
the  zero  order  problem  given  in  chapter  VII  B.  This  dependence  can¬ 
not  be  neglected  and  consequently  we  have  to  use  for  these  terms  new 
developments  of  the  form 

*‘V'  -  a,,  J,  (^r/R) 

“63)  rf'%  -  I  Je  (/9„  r/R) 


relation  between  the  corresponding  coefficients  and  n^vand  T 

/**,v 
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(165)  -  —  JjOh, T]T  ^  r/R)  JC r/R)rdr 

Now  we  introduce  equations  (160),  (162),  (164)  and  (165)  into  equations 
(157)  and  (158)  and  collect  all  terms  with  the  same  value  Jq  (J5y  r/R). 

Then  the  coefficients  of  these  terms  must  all  be  zero  due  to  the  fact 
that  the  relations  (157)  and  (158)  hold  for  any  value  of  r  in  the  range 
o  &  r  £  R.  This  produces  an  infinite  system  of  coupled  linear  equa¬ 
tions  for  the  coefficients  n  vand  T  .  Since  we  have  as  many  homo- 
gaieous  equations  as  we  have  variables,  the  determinant  of  the  system 
must  vanish.  This  secular  determinant  is 


MG 


■‘SW#  M'-iV- 


87 


CO  f  "  Cc  fl0#  j  -  c0  B,c  J 

Si 

O 

• 

-  C,  B01  ;  to+ysfi -C.Bn;  •  • 

• 

'  *  CiAsi  •  *  • 

••(*?-)*  i  0 

0 . 

0  i  -‘M'1/ 0  ” 

• 

• 

R 

S«  •/•‘"’-J./ r/R)rdr 

0 

r  .  2 

R 

At*  "  f  “Mo  ft,  'AM&  r/#^rdr 

s 

The  secular  determinant  defines  the  eigenvalues  of  our  system.  Intro¬ 
ducing  one  of  these  eigenvalues  in  the  equations  allows  to  calculate  the 
coefficients  of  the  corresponding  eigensolution. 

Our  stability  problem  is  therewith  reduced  to  the  solution  of  the  secular 
determinant. 

There  are  two  reasons  why  this  solution  is  extremely  difficult.  One  is 
that  we  are  dealing  with  an  infinite  determinant,  the  other  that  the  coef¬ 
ficients  occuring  in  this  determinant  are  complicated  integral  expressions. 
Approximations  will  therefore  be  unavoidable. 

We  suggest  to  simplify  the  problem  by  omitting  the  terms  -V  >  <8 
in  the  relations  (164,165).  This  seems  possible  since  the  numerical 
value  of  the  integral  expressions  decreases  rapidly  with  increasing 
index  number  v  **  •  .  This  is  particularly  true  for  small  values 

of  s.  It  is  favourable  that  our  investigation  is  just  interested  in 
these  eigensolutions  with  a  small  number  of  modes  (small  s)  which 
decide  the  onset  of  the  instability. 
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Under  these  circumstances  the  determinant  (166)  after  suitable  rear¬ 
rangement  can  be  written  as  a  product  of  determinants  in  the  form 


-C,  ,  Cv  4, 


(167)  ^  •  H 


Each  of  these  factors  provides  a  quadratic  equation  for  two  of  the 
possible  eigenvdues  .  From  chapter  III  we  know  that  a  necessary 
and  sufficient  condition  for  stability  is  that  the  coefficients  in  this 
relation  have  the  same  sign. 

This  provides  for  the  zero  order  mode  the  following  instability  criterion 

Ota)  X  n*  3»a  r._  ills.)] 
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The  higher  order  modes  become  unstable  if  the  condition 


is  fulfilled. 


Here  we  have  used  the  abbreviations 
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An  estimation  of  the  quantities  entering  our  criterion  shows'that  the 
onset  of  the  instability  is  determined  by  the  zero  order  mode. 

Since  our  calculations  do  not  exclude  the  case  of  the  normal  glow  column, 
let  us  investigate  whether  our  criterion  (168)  supports  our  expectation 
that  the  normal  glow  column  should  be  always  stable. 

In  the  case  of  S  c  h o  t  tky "  s  theory  di 
discharge  volume.  Consequently  we  find 
the  eigenvalue  «6^  /D  ^  =  (2,4/R)2 
arrive  at 

and  should  ’conclude  -  contrary  to  our  expectations  -  that  the  normal 
glow  column  is  always  unstable.  This  conclusion  is  of  course  incorrect. 
Since  in  the  normal  Schottky  theory  the  first  term  in  equation  (168) 
is  small  of  higher  order  we  must  include  in  our  stability  consideration 
higher  order  contributions  to  the  stationary  eigenvalues 

Let  us  therefore  study  higher  approximations  of  the  stationary  solution 
by  using  the  development 
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Using 

from  chapter  VII  A  we 
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Introducing  this  into  equations  (55,63)  we  find 
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Here  is  the  zero  order  of  the  eigenvalue  (<*(oo))  =(a>/«)*0«f 

s"'i  8  the  correction  term  of  the  eigenvalue.  oC(°D 

is  the  perturbing 

function  which  can  be  found  from  equation  (30)  using  the  relation 
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From  well  known  perturbation  theory  it  follows  therefore  for  the 
first  order  correction  of  the  eigenvalue  oC^00) 
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Introducing  in  our  criterion  (23)  gives  after  some 

simple  rearrangements  the  condition  for  instability  1  »  1.  Clearly 
this  condition  cannot  be  fulfilled  and  as  expected  the  Schottky  column 
is  always  stable.  If  we  carry  the  perturbation  theory  on  to  the  second 
order  this  fact  is  emphasized. 


The  situation  is  quite  different  for  the  thermal  glow  column.  We  refer 
to  the  results  calculated  in  chapter  VII  B.  Introducing  the  eigenvalue 
found  in  equation  (81  )  our  criterion  determining  the  stability  reads  now 


<*&**&*■*>, 
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and  we  see  that  instability  of  the  thermally  constricted  Schottky 
column  can  well  occur. 


IX  INSTABILITY  OF  THE  CATHODIC  PART 
OF  THE  GLOW  DISCHARGE. 


In  the  following  we  discuss  a  special  type  of  instability  of  the  cacho- 
dic  part  of  the  glow  discharge:  The  "  8*  -  <f  instability". 

Let  us  recall  the  phenomenology  and  the  mechanism  of  the  cathodic 
part  of  the  glow  discharge. 

The  cathodic  part  of  the  glow  discharge  is  that  region  which  is  in¬ 
fluenced  by  the  cathode.  If  we  carry  out  an  infinitesimal  displace¬ 
ment  of  the  cathode  keeping  the  current  constant  then  all  regions 
where  we  find  changes  in  the  discharge  parameters  belong  to  the 
cathodic  part. 

The  so  defined  cathodic  part  of  the  glow  discharge  is  subdivided  in  a 
number  of  phenomenological  zones  and  a  number  of  model  regions. 
The  two  subdivisions  are  not  identical. 

The  phenomenological  regions: 

1)  The  Aston  dark  space  is  in  general  small  and  may  be  absent  with 
high  cathode  falls.  It  lies  immediately  in  front  of  the  cathode  and 
is  characterized  by  the  fact  that  the  electrons  deliberated  from  the 
cathode  do  not  have  sufficient  energy  to  cause  notable  excitation. 

2)  The  first  cathode  layer  is  a  small  layer  of  faint  light  which  is  dif¬ 
fuse  in  the  direction  towards  the  anode.  This  is  a  zone  in  which 
the  electrons  have  just  reached  energies  corresponding  to  the 
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first  excitation  levels  of  the  gas.  Its  decreasing  visibility  to¬ 
wards  the  anode  is  due  to  the  acceleration  of  the  electrons  and 
the  corresponding  decrease  in  their  excitation  probability. 

3)  HittorPs  or  Cfrook^s  dark  space  appears  dark  in  contrast  to  the 

neighbouring  regions  of  the  first  cathode  layer  and  the  negative 
glow.  When  the  electrons  enter  the  Hittorf  or  Crook  dark  space 
they  have  already  energies  above  the  excitation  energies  of  the 
gas  and  accummulate  further  energy  in  the  electric  field.  Due  to 
this  influence  they  finally  surpass  the  state  of  drift  motion  and 
become  runaways  .  The  beam  character  of  the  electrons 

in  the  cathode  fall  region  has  been  proved  by  a  number  of  experi¬ 
ments.  One  of  the  most  obvious  proofs  is  provided  by  an  obstacle 
in  the  Crook  dark  space.  Under  these  circumstances  the  part  of 
the  negative  glow  behind  the  obstacle  disappears  and  also  all  the 
corresponding  cathode  parts  in  the  area  of  the  obstacle  are  des¬ 
troyed. 

4)  The  negative  glow  is  the  most  luminous  region  of  the  whole  cathode 
zone  with  a  sharp  boundary  on  the  cathode  side.  This  boundary  is 
called  the  glow  edge.  There  is  a  continous  decrease  in  brightness 
towards  the  anode.  The  geometry  of  the  negative  glow  is  dependent 
only  on  the  geometrical  form  of  the  cathode  and  independent  of  po¬ 
sition  and  form  of  the  anode.  This  is  due  to  the  fact  that  the  elec¬ 
tron  beams  are  practically  always  perpendicular  to  the  cathode  sur¬ 
face. 

5)  Faraday  dark  space  is  a  region  of  low  luminosity  which  is  limited 
by  the  anode  itself  or  by  the  head  of  the  positive  column.  The  length 
of  the  Faraday  dark  space  increases  with  the  radius  of  the  discharge 
vessel  and  is  in  general  larger  than  the  extension  of  the  negative  glow. 
In  the  Faraday  dark  space  the  charged  particles  diffuse  towards  the 
anode  and  towards  the  wall.  There  is  virtually  no  electric  field  and 

no  particle  production  in  this  region. 
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The  model  regions  describing  the  situation  in  the  cathodic  part  of  the 
glow  discharge  are  the  cathode  fall  region  (CFR),  the  nega¬ 
tive  glow  region  (NGR)  and  the  F  a  r  a  d  a  y  dark  space  re¬ 
gion  (FDR).  We  have  already  discussed  the  properties  of  these  model 
regions  in  our  description  of  the  s-mode  of  the  cathodic  part  of  the 
glow  discharge  in  chapter  VII  I  . 

It  is  impracticable  to  treat  the  complicated  problem  of  instability  of 
the  cathode  part  of  the  glow  taking  into  account  all  the  details  and  ex¬ 
tensions  given  in  FTR  II.  Rather  we  will  have  to  restrict  ourselves 
here  from  the  very  beginning  to  tne  most  simple  description  of  the 
cathodic  part  of  the  glow  discharge,  which  is  possible  without  affecting 
the  principal  features  of  the  instability  criterion. 

It  is  sufficient  to  consider  the  CFR.  The  conditions  in  the  negative 
glow  enter  these  calculations  via  the  boundary  condition 

}f(d)  •$}.(<*)  0?9) 

A  complete  theory  would  have  to  calculate  the  coefficient  &  from 
the  theory  of  the  NGR.  This  was  not  achieved  in  the  case  of  the  sta¬ 
tionary  solution  and  is  of  course  even  less  possible  in  the  nonstationary 
case.  However  this  lack  of  knowledge  will  not  affect  our  following  cal¬ 
culations  since  we  allow  quite  generally  for  a  variation  of  S  with  the 
electron  current  density  j_(d)  . 


We  write  down  the  nonstationary  equations  governing  the  CFR  in  the 
usual  one-dimensional  approximation.  We  have  the  particle  conserva¬ 
tion  laws 


(180) 


(181)  -  O 

where  j_,  j  describes  the  current  density  and  n_,  n+  the  average 
particle  number  density  for  the  electrons  and  ions  respectively. 

The  coefficient  06  designates  the  production  of  particles  per  cm  and 
electrons  in  the  CFR.  This  is  in  general  a  very  complicated  quantity 
since  the  electron  in  the  CFR  have  partly  runaway  character.  Their 
velocity  depends  on  their  origin  in  the  CFR  and  consequently  their 
velocity  distribution  is  upto  now  unknown.  The  problem  has  been  dis- 
cussed  in  FTR  II.  In  general  tb  must  be  considered  as  a  function  of 
the  electric  field  E{z)  and  the  time  t 


<*'  I'M! 


However  we  will  study  a  variation  of  our  system  with  a  characteristic 
time  large  in  comparison  to  the  relaxation  time  of  the  electron  motion 
and  production.  Under  these  circumstances  the  motion  of  the  electrons 
and  their  production  is  defined  by  the  momentary  electric  field  E(z) 


in  the  CFR  and  we  have 


o6*(E) 


and  for  the  momentum  balances 


-  at.  (e)  n. 
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r-f  1  V*6#r-MVv~  H-.*a»o»i,.i«»>* 


95 

where  the  coefficient*  M.  end  ere  both  positive  quentitie*.  The 
geometry  of  our  coordinete  system  is  shown  in  Fig.  22 

For  the  ions  it  may  be  sefe  to  assume 
mobility  limited  motion.  For  the  elec¬ 
trons  this  assumption  is  rather  doubt¬ 
ful.  Fortunately  we  do  not  need  to  spe 
cify  the  quantities  >€„  86f  . 

We  have  further  the  Poisson  equa¬ 
tion 

-  tf(n4-nj)  (186) 

Equations  (180,181)  and  (183-186)  define  the  situationin  the  CFR  together 
with  the  boundary  conditions  at  the  cathode  surface  and  at  the  anode  sur¬ 
face. 

At  the  NGR  boundary  we  have  the  condition  (179).  We  recall  that  the 
coefficient  &  in  general  must  be  considered  depending  on  the  current 
j_(d).  Also  at  the  NGR  boundary  we  assume  for  the  electric  field 

B(d)  -  o  <187) 

At  the  cathode  surface  the  electrons  deliberated  from  the  cathode  are 
due  to  ion  bombardment  of  the  cathode  surface.  Therefore  we  have 
the  condition 


e 

■  f 

T  * 

rig.  22 1  Coordinate  system  and  symbols  used 
in  tbs  cathodic  part  of  the  glow. 


j,  Jo)  •  ?  (188) 

where  the  coefficient  if  -  for  the  same  reason  as  S  •  may  be  still 
a  function  of  the  current  density  j+(o). 
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We  auggeat  to  study  the  6  instability  defined  by  the  following  con- 
cept. 

We  have  a  given  atationary  aolution  of  the  model  regions  of  the  catho¬ 
dic  part  of  the  glow  discharge  as  described  in  FTR  II.  It  is  assumed 
that  radial  boundary  effects  are  negligible,  so  that  the  whole  situation 
can  be  treated  one-dimensional.  Now  we  superimpose  a  momentary 
disturbance  of  the  parameters  describing  the  cathodic  part.  This 
disturbance  is, limited  to  a  small  fraction  of  the  whole  cathode  area. 

The  lateral  extension  of  this  disturbance  is  subject  to  two  counteracting 
conditions.  On  the  one  hand  the  lateral  extension  is  still  large  enough 
to  allow  the  one-dimensional  description.  On  the  other  hand  the  lateral 
extension  of  the  channel  has  to  be  small  enough  so  that  most  of  the 
field  lines  originating  from  the  space  charge  caused  by  the  disturbance 
leave  the  discharge  channel  in  lateral  direction  and  do  therefore  not 
effect  the  coefficients  H  and  0 6*  .  There  is  a  wide  range  of  dis¬ 
turbances  for  which  these  conditions  are  fulfilled  simultaneously. 


To  the  so  defined  model  we  apply  perturbation  theory  and  have 


n  ■  n  *-  n 

m 

• 

0 

(189) 
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Here  the  index  (o)  characterizes  the  stationary  solution,  the  index  (1) 
the  nonstationary  disturbance. 


Recalling  our  assumptions  introduced  above  the  first  order  components 
of  the  current  densities  are  defined  by  the  equations 


(190) 


dt  3z 


♦ 
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and 


iU.'/Kj  Zlif’/a)  ..  M  . 

it  »z  *  lt./«l  m  o  m 

The  boundary  conditions  (179)  and  (188)  read  in  the  first  order  approxima¬ 
tion,  -  remembering  the  dependence  of  6  on  j_(d)  and  T  on 

j+  (°) 

t,'ld)  -  i ^ ld)  ,  <1,2) 


*  ir ,  r  •  ru>t  j%M  <•>  (,,3) 

We  look  for  the  eigensolutions  of  this  problem  by  introducing 

f’u.t)  -  j'*W*  <>’«> 

into  equations  (190)  and  (191).  The  result  is 

j 
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In  these  equations  the  symbol  Sg  is  identical  with  the  sign  of  the 
current  .  Equation  (195)  may  be  at  once  integrated  resulting  in 


/l*1 &)«<*} 
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where  we  have  now  introduced  a  new  production  coefficient  of  the  form 


(198) 


06  m  oC  Sg 


(A) 


This  fact  should  be  carefully  observed  in  case  j_  goes  through  zero 

The  method  of  the  variation  of  the  constant  gives  for  the  second  equa¬ 
tion 

.<»>  •  * 
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(199) 
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Applying  the  boundary  conditions  (193)  and  (192)  we  find 


(200) 


and 


(201) 


i^(o)  -  jl»(o)/r 


f  mxp  j- 
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Since  Af.  is  very  much  larger  than  we  may  reduce  this 

equation  to 

(202)  <f  {■-  J  ^ )  dx J  -  [  £  ♦  j  5  {  J  (**  )  ftt} J 

1  •  #o 


We  further  introduce  the  new  variable 


04 


U 

which  gives 


0 


u  (d)*  to 


(203) 


Within  the  frame  of  our  approximations  it  is  quite  safe  here  to  consider 
the  term  Z  /  (t i  *  in  the  integral  (204)  as  constant  in  compa¬ 

rison  to  the  exponential.  We  replace  it  by  the  average  value 


r-  d  -  d- 

j  Zd*j  ^  (Z+ ^)dz  -  j Zdz/uJ  (205) 
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Our  cathodic  part  of  the  glow  discharge  can  be  considered  stable 
then  and  only  then  if  the  real  parts  of  the  complex  numbers  X^ 
have  all  negative  values.  With  other  words  the  relation  must  hold 

d 

(208)  a,;  <  j*o/z 

0 

Let  us  now  -  following  Nyquist  -  consider  the  function 


(209) 


d 

JV  -  *  foZc/z/u)Jexp(-u>)  - 


d 

jtCdz/u) 


which  describes  a  transformation  from  the  to  -  complex  plane  into 
the  W- complex  plane. 


We  decide  on  the  question  of  stability  by  using  the  relations  (7)  and 
(8)  of  chapter  III.  According  to  equation  (209)  we  have  to  choose  our 
integral  contour  for  equation  (7)  as  shown  in  Fig.  23.  (remember 
•  that  is  a  negative  quantity  due  to  the  sign  j_  ). 


To  see  under  what  circumstances  the  path  in  the  W-plane  encircles 
the  origin  let  us  construct  the  W  number  vector  by  adding  it  up 
from  the  various  contributions  of  relation  (209). 

There  is  first  the  constant  contribution  -1/  y-  .  (See  Fig.  23). 

The  second  contribution 

(210)  •  joidz/u)  »  JoC  dz/(x  *  iy) 

•  0 

is  in  comparison  to  the  exponential  contribution  only  little  varying. 
Its  absolute  value  is  also  small  in  comparison  to  the  exponential 
contribution.  It  is  given  in  Fig.  23  by  the  vector  .  Finally 

there  is  the  third  contribution 


r  gf* 


1 
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_  c { 

^  S  *  JoC  difu>  J  •up  (-at)  »  f  r  (211) 

consisting  of  two  parts  which  are  represented  in  Fig.  23  by  the 
vectors  and  .  Jr  rotates  around  the  point  -Vff, 

in  the  W-plane  with  a  frequency  defined  by  the  imaginary  part  of  u)  . 

A  detailed  discussion  of  the  vector  diagram  of  Fig.  2  3  indicates  that 
the  largest  value  of  the  real  part  of  W  occurs  for 

d 

y  •  o  and  x  ■  -  J  jot  J  dt  (212) 

0 

Therefore  encircling  of  the  W-origin  can  be  expected  only  if  for  the 
values  of  equation  (212)  W  has  a  positive  real  part  (Wr  >  o).  It 
follows 

.  d 

(j  +  6  )  J  l<*l  dzj  •<.  7*4  (213) 

o 

for  the  stability  of  our  system. 

Let  us  remember  that  the  stationarity  condition  has  the  form 

d 

(/*6)e*p|  Jk|c/*  }-  1  +  -?  (214) 

0 

Then  we  recognize  that  ^  -  6  instability  will  occur  if  the  relation 


f  (ru)(f  »i)  y  1 
f  (f-i)  (.6  *  i ) 


holds. 


102 


If  either  <J  or  JT  or  both 
of  these  quantities  are  larger 
than  S  or  Jf  ,  then  in¬ 
stability  will  occur.  But  insta¬ 
bility  can  even  occur  if  one  of  - 

these  quantities  is  sufficiently  fat 
larger  than  the  undisturbed 
value  and  the  other  one  is  small 
but  in  effect  not  overcoming  the 
influence  of  the  first  quantity. 


Pig.  23s  Ny quiet  mapping  for  the  dispersion 
relation  of  the  instability. 


These  results  are  quite  visible.  If  an  increase  in  the  current  density 
of  the  electrons  results  in  an  ion  component  from  the  negative  glow 
larger  than  that  necessary  for  the  stationary  state  we  will  expect  an 
avalanche  process  causing  instability.  The  same  holds  if  an  increase 
in  the  ion  current  density  produces  an  electron  cathode  emission  too 
large  for  the  stationary  state. 

It  should  be  noted  particularly  that  a  spontaneous  increase  in  the  cur¬ 
rent  density  of  the  CFR  itself  does  not  produce  instability.  The  insta¬ 
bility  occurs  only  if  the  current  density  increase  goes  along  with  a 
change  in  the  coefficients  and  S  ! 

It  is  very  difficult  to  relate  these  results  to  experimental  observations. 
There  are  quite  a  number  of  experiments  (e.  g.  91,92,93,58,18,19,25.27,28, 
95,  96,  8, 10,  38  )  which  show  cathodic  glow  instability  and  point  directly 
to  the  importance  of  small  size  surface  impurities.  The  experiments 
would  allow  an  explanation  with  our  2T  -  6  instability  but  they  do  not 
exclude  other  possibilities. 
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X  STABILITY  PRINCIPLES  IN  GASEOUS  ELECTRONICS. 

In  chapter  III  we  sketched  a  general  method  to  derive  variational  prin¬ 
ciples  for  stability.  This  method  requires  a  certain  set  of  equations 
which  describes  the  whole  phenomenon  under  consideration.  As  an 

example  we  refer  to  the  calculation  of  the  magneto  hydrodynamic 
q 

principle 

The  application  of  this  general  method  to  problems  of  gaseous  elec¬ 
tronics  confronts  us  with  the  difficulty  that  there  is  in  gaseous  elec¬ 
tronics  no  unique  set  of  differential  equations  describing  the  whole 
field.  Rather  in  any  gas  discharge  we  have  a  certain  number  of  model 
regions  and  each  of  them  is  governed  by  different  sets  of  equations. 

In  contrast  to  the  situation  in  magneto  dydrodynamics  we  can  there¬ 
fore  not  expect  to  find  one  principle  governing  all  gas  discharges. 

In  the  best  case  we  can  hope  to  have  different  principles  for  different 
discharges  and  discharge  parts. 

The  derivation  of  such  variational  principles  is  a  very  involved  problem 

q 

as  can  be  seen  for  instance  from  the  investigation  7  .  It  is  there¬ 

fore  beyond  the  scope  of  this  investigation. 

We  will  however  discuss  here  two  important  problems.  The  one  is  the 
question  question  whether  the  well  known  "principle  of  minimum 
entropy  production”  of  irreversible  thermdynamics  may  be  used 
to  decide  on  the  stability  in  gaseous  electronics.  The  second  is  the  dis¬ 
cussion  of  "S  t  e  e  n  b  e  c  k' s  minimum  principle"  which  has  been 
widely  applied  in  gaseous  electronics  to  decide  on  stable  modes. 

A  THE  PRINCIPLE  OF  MINIMUM  ENTROPY  PRODUCTION. 

The  "principle  of  minimum  entropy  production"  can  be  applied  if  the 
following  conditions  are  met: 
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1.)  Gibb's  fundamental  equation 


(216)  TdS  -  dU  f  Pd V 


must  hold. 


2  a,  dMi 


2.)  The  linear  phenomenological  relations  between  the  thermodynamical 
fluxes  and  forces  should  be  valid. 


3. )  The  Onsager  Casimir  reciprocity  relations  should  be  appli¬ 
cable. 


4. )  The  phenomenological  coefficients  should  be  time -independent. 

5.  )  The  system  should  be  in  mechanical  equilibrium  or  at  least  in  a 
stationary  mechanical  state  according  to 


(217) 


bt 
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6. )  At  all  parts  of  the  surface  of  the  system  one  of  the  two  following 
conditions  must  hold.  Either  the  temperature  and  the  chemical  poten¬ 
tials  should  be  time -independent  or  the  temperature  should  be  time- 
independent  and  the  normal  component  of  all  mass  current  should  dis¬ 
appear. 


These  conditions  imply  that  each  component  of  the  system  is  close  to 
equilibrium.  The  principle  of  minimum  entropy  production  states: 

The  stationary  states  of  a  system  subject  to  the  above  conditions  are 
characterized  by  a  minimum  of  entropy  production.  These  states  are 
also  stable 
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If  we  write  down  the  equations  of  irreversible  thermodynamics  for  a 
plasma  satisfying  the  above  conditions  and  apply  the  principle  of  mi¬ 
nimum  entropy  production  then  the  Euler  equations  of  the  varia¬ 
tional  problem  are  identical  with  the  transport  equations.  Since  we 
are  on  the  other  hand  close  to  equilibrium  the  Saha  equation  is 
valid.  The  application  of  the  energy  transport  equation  and  the  Saha 
equation  is  the  basic  scheme  of  the  "arc  approximation"  described 
in  chapter  VI  B. 

Therefore  an  important  conclusion  can  be  drawn  from  the  principle 
of  minimum  entropy  production: 

Any  discharge  part  which  can  be  described  by  the  arc  approximation 
and  has  boundary  conditions  which  meet  the  conditions  (6)  and  (7)  is 
always  a  stable  one. 

On  the  other  hand,  any  discharge  part  which  is  not  governed  by  the 
laws  of  the  arc  approximation  cannot  be  subject  to  the  application 
of  the  principle  of  minimum  entropy  production. 

Therefore  the  wall  stabilized  arc  column  treated  in  chapter  VII  E 
is  always  a  stable  mode.  By  the  same  token  all  other  modes  described 
in  chapter  VII  cannot  be  subject  to  the  application  of  the  principle  of 
minimum  entropy  production. 

However  it  is  still  possible  to  apply  the  principle  to  that  part  of  a  dis¬ 
charge  mode  described  by  the  arc  approximation  -  provided  that  the 
interface  limiting  this  region  is  not  varied. 

With  this  restriction  the  thermal  nucleus  of  the  wall  detached  arc 
column  shows  the  same  stability  as  the  wall  stabilized  arc. 
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If  we  allow  for  variations  of  the  interface  between  the  two  model 
regions  then  no  conclusion  on  the  stability  can  be  drawn  from  the 
application  of  the  minimum  entropy  production  principle  since  all 
the  calculations  carried  out  for  this  principle  underlie  the  basic 
assumption  that  we  have  boundary  conditions  fixed  in  space.  The 
latter  problem  of  the  stability  of  the  wall  detached  arc  column  with 
respect  to  variations  of  the  diameter  of  the  thermal  nucleus  leads 
us  over  to  the  problem  of  the  Steenbeck  minimum  principle. 


B  Steenbeck"  a  MINIMUM  PRINCIPLE. 

Following  an  investigation  of  Compton  and  Morse, 

Steenbeck  has  claimed  the  following  postulate: 

In  a  real  gas  discharge  of  constant  total  discharge  current  the  energy 
production  and  with  that  the  voltage  applied  to  the  electrodes  is  a 
minimum. 

If  one  identifies  a  "real  discharge"  with  a  discharge  which  is  stable 
then  Steenbeck's  minimum  principle  is  a  variational  principle 
for  stability. 

Steenbeck  's  minimum  principle  could  be  varified  in  a  number  of 
experiments,  e.  g.  20,83,99  .  The  theoretical  discussion  has  in¬ 

vestigated  whether  Steenbeck's  empirical  postulate  can  be  justified. 

There  are  essentially  two  questions. 

1)  Is  Steenbeck's  minimum  principle  a  true  variational  principle 
for  stability? 

2)  If  it  is  a  true  variational  principle  what  are  the  conditions  for  the 
application  of  this  principle? 
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70 

W  e  i  z  e  1  and  R  o  m  p  e  investigated  the  question  whether  in  the 

range  of  the  arc  approximation  the  Elenbaas-Heller  equation  is 
Euler's  equation  of  S  t  e  e  nb  e  c  k'  s  minimum  principle.  They  come 
to  the  conclusion  that  in  general  this  is  not  the  case. 

65 

Peters  starts  from  the  minimum  entropy  production  principle 

of  irreversible  thermodynamics  and  tries  to  deduce  Steenbeck's 
postulate  from  this  principle.  He  finds  Steenbeck's  principle. 
However  K  r  a u  s  e  andSteenbeck  ^  pointed  out  that  these  cal¬ 
culations  include  the  inaccurate  restriction  to  variations  into  stationary 
states. 

According  to  K  r  au  8  e  andSteenbeck  such  variations  into  sta¬ 
tionary  states  are  possible  if  one  uses  the  "canal  model".  Therefore 
for  this  canal  model  St  e  e  nb  e  c  k' s  minimum  principle  follows  from 
the  principle  of  minimum  entropy  production  and  for  this  special  case 
is  therefore  manifested  as  a  true  variational  principle. 

20 

In  consideration  of  the  canal  model  other  investigators 
came  to  the  conclusion  that  Steenbeck's  principle  is  a  mean  to 
determine  the  additional  parameter  of  the  discharge  diameter  intro¬ 
duced  by  the  canal  model. 

The  most  refined  discussion  of  Steenbeck's  principle  has  been  given 

84 

recently  by  Ste  enbe  ck  himself  .  Steenbeck  states  that 

Elenbaas-Heller  's  equation  is  only  an  approximation.  A  more 
exact  description  including  additional  effects  would  increase  the  num¬ 
ber  of  degrees  of  freedom.  The  corresponding  boundary  conditions  for 
these  variables  may  be  ineffective.  Then  Steenbeck  says  there  is 
not  only  space  for  the  application  of  a  minimum  principle  but  there  is 
even  the  necessity  of  the  application  of  such  a  principle  to  select  the 
actual  discharge  mode. 

Our  interpretation  of  the  Steenbeck  principle  may  be  described  as 
follows: 
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Let  ua  consider  quite  generally  any  kind  of  a  discharge  volume  des¬ 
cribed  by  a  certain  set  of  differential  equations.  Let  us  assume  further 
that  we  have  a  sufficient  number  of  effective  boundary  conditions  for 
all  variables  occuring  in  these  equations. 

Then  there  is  one  and  only  one  stationary  solution  for  these  differential 
equations  with  the  given  boundary  conditions.  If  there  is  at  all  a  stable 
solution  it  must  be  this  stationary  solution. 

Now  any  physical  system  is  subject  to  statistical  fluctuations.  If  the 
influence  of  one  or  more  of  the  boundary  conditions  is  decreased  then 
we  reach  a  point  where  the  influence  of  these  boundary  conditions  comes 
below  the  level  of  the  average  statistical  fluctuations.  In  this  case  the 
corresponding  boundary  condition  becomes  ineffective  (see  chapter  VI  C). 

Due  to  the  lack  of  this  boundary  condition  we  do  not  have  only  one  solu¬ 
tion  but  a  whole  manifoldness  of  stationary  solutions  of  our  system.  In 
a  physical  experiment  we  will  observe  that  one  of  all  these  stationary 
solutions  which  is  stable.  Under  these  circumstances  we  therefore  al¬ 
ways  need  a  stability  criterion  which  selects  the  solution  observed  in 
the  experiment.  The  kind  of  this  criterion  depends  on  the  physical  si¬ 
tuation. 

With  respect  toSteenbeck's  minimum  principle  we  now  consider  a 
discharge  consisting  of  an  arc  model  region  and  a  glow  model  region 
separated  by  an  interface  (e.  g.  wall  detached  arc  column,  cathodic 
part  of  the  arc  discharge). 

As  was  explained  in  the  chapters  VI  C,  VII  F,  VII  G  the  boundary  con¬ 
ditions  for  this  arrangement  may  become  ineffective  if  the  walls  are 
removed  too  far.  Then  the  shape  of  the  discharge  is  not  uniquely  deter¬ 
mined  anymore.  In  our  model  representation  this  is  reflected  in  the 
variability  of  the  position  of  the  interface.  For  the  case  of  the  wall 
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detached  arc  column  the  poeition  of  the  interface  is  described  by  the 
diameter  of  the  arc  nucleus. 

It  is  our  opinion  that  Steenbeck's  minimum  principle  selects 
that  mode  which  is  stable  with  respect  to  disturbances  in  the  position 
of  the  interface.  With  respect  to  all  other  disturbances  the  thermal 
nucleus  is  stable  according  to  the  principle  of  minimum  entropy 
production.  Therefore  -  for  our  special  arrangement  of  an  arc  model 
region  and  a  glow  model  region  -  Steenbeck's  minimum  principle 
is  a  sufficient  criterion. 


"Jr***- 
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